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ABSTRACT 


TVHIS, a linear time-varying helicopter 
simulation program, is described in this report. 

The program is designed as a realistic yet effi- 
cient helicopter simulation. It is based on a 
linear time-varying helicopter model which in- 
cludes rotor, actuator, and sensor models, as 
well as a simulation of flight computer logic. 

TVHIS can generate a mean trajectory simulation 
along a nominal trajectory, or propagate covari- 
ance of helicopter states, including rigid-body, 
turbulence,- control command, controller states, 
and rigid-body state estimates. 

Rotor and actuator models include dynamic 
states with actuator displacement and rate limits. 
Flight computer logic includes state estimation, 
guidance, and digital-adaptive control laws. Auto- 
matic guidance may be specified to cause the heli- 
copter to follow a prestored guidance trajectory, 
or pilot input may be applied as perturbations on 
nominal guidance commands. State estimation and 
guidance algorithms may be by-passed to allow the 
testing of control algorithms alone. 
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1 . 


TVHIS OVERVIEW 


TVHIS is designed to simulate the flight path of a 
helicopter in the neighborhood of a reference flight trajec- 
tory. Either a mean trajectory simulation, or a covariance 
calculation may be performed. Multiple cases may be run, 
combining different segments of a nominal trajectory, dif- 
ferent types of control and guidance for a given trajectory 
segment, or covariance calculations along different segments 
of the nominal trajectory. 

TVHIS is functionally divided into 6 sections: 

(1) Main Executive Controlling Routine 

(2) Data Base Modification and Echo 

(3) Nominal Trajectory Construction 

(4) Propagation Between Way Points 

(a) State Equation Integration 

(b) Covariance Equation Integration 

(5) Flight Computer Logic 

(6) Time History Summary 

(a) State Time History 

(b) Covariance Time History 

When TVHIS is first loaded, the TVHIS data base, con- 
sisting of general variables and matrices used for control 
and flight computer variables, is initialized to default values 
specified by two BLOCK DATA subprograms. One subprogram speci- 
fies defaults for general variables, and a second subprogram 
specifies defaults for flight computer control settings. 
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R— 31819 



Figure 1-1 TVHIS Data Flow 


During execution, TVHIS reads input data from a 
card-image input file that describes the specific run being 
performed (Fig. 1-1). Aerodynamic coefficients are read 
from a permanent disk file, which must be established using 
the procedure described in Appendix B, before the execution 
of TVHIS. 

1.1 MAIN EXECUTIVE ROUTINE 

The main routine (Fig. 1.1-1) reads namelist 
$STACK, specifying NCASES , the number of cases required in 
the run. For each of the NCASES steps, the following will 
be performed: 
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Figure 1.1-1 TVHIS Main Executive Block Diagram 
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(1) Data Base Modification and Echo (RDDAT) 

(2) Nominal Trajectory Construction (MISION) 

(3) Wind and Sensor Initialization During 
First Flight Segment for a Mean Trajec- 
tory Simulation ( SENINT, WNDINT) 

(4) Propagation Between Nominal Waypoints 

(a) State Equation Integration Between 
Waypoints for Mean Trajectory Simu- 
lation (RKGS) 

(b) Covariance Equation Integration 
Between Waypoints for Covariance 
Calculator (INTEG) 

(5) Final Time History Print and/or Plot 

(a) Mean Trajectory Simulation (MSUMRY) 

(b) Covariance Calculation (PSUMRY) 


1.2 DATA BASE MODIFICATION AND ECHO 

The data input subroutine RDDAT (Fig. 1.2-1) reads 
modifications to the TVHIS data base via NAMELIST and issues 
a report on the initial data base at the beginning of each 
case. Inputs read by RDDAT are described in Section 2.1.1; 
outputs are described in Section 2.2.1. In the case of a co- 
variance calculation, RDDAT calls COVPT to produce a report 
of the initial values of covariance matrices. 


1.3 NOMINAL TRAJECTORY CONSTRUCTION 

Subroutine MISION (Fig. 1.3-1) calculates system 
matrices F, G, HT, and nomin al state XVNTRU at two waypoints 
of the nominal trajectory which bound the current flight seg- 
ment . 
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RDDAT 


DATA BASE 
MODIFICATION 
AND ECHO 


COVPT 


PRINT 

NORMALIZED 

COVARIANCE 

MATRICES 


Figure 1.2-1 RDDAT Block Diagram 
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Figure 1.3-1 MISION Block Diagram 
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FLTINT interpolates values from aerodynamic coef- 
ficient data to determine system matrices. Aerodynamic 
coefficients are read from a permanent disk file on FORTRAN 
logical unit 1. If required, DYNTRM is called for circular 
and accelerating trajectories to calculate dynamic trim 
corrections and dynamic trim Euler angles. 


1.4 STATE EQUATION INTEGRATION 

RKGS (Fig. 1.4-1) was adapted from the IBM Scientific 
Subroutine Package (SSP) (Ref. 1). RKGS numerically solves a 
system of first-order differential equations using the fourth 
order Runge-Kutta Method with modifications due to Gill. The 
integration step-size is automatically adjusted to produce the 
desired degree of accuracy. 

Subroutine MCON evaluates the derivatives of con- 
tinuous states. MDIG records values of continuous states, 
propagates discrete states, executes sensor transformations, 
and calls the flight computer logic. MDIG calls MSUMRY if 
matrices used for intermediate time-history storage are full. 


1.5 COVARIANCE EQUATION PROPAGATION 

Subroutine INTEG (Fig. 1.5-1) numerically integrates 
a vector representing a covariance matrix by rectangular or 
trapezoidal integration techniques. PCON is called to cal- 
culate the derivative of covariance of continuous states. 

PDIG updates the nominal trajectory and calls subroutine 
PCMPTR to propagate covariance of discrete states from the 
flight computer. PDIG calls PSUMRY if intermediate time- 
history storage matrices are full. 
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Figure 1.4-1 RKGS Block Diagram 
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Figure 1.5-1 


INTEG Block Diagram 


1.6 FLIGHT COMPUTER LOGIC 

Subroutine CMPUTR (Fig. 1.6-1) is designed to follow 
the helicopter flight computer logic described in References 
2, 3 and 4. 
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Figure 1.6-1 


CMPUTR Block Diagram 
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The basic functions of CMPUTR, and the subroutines 
performing these functions are: 


(1) Reset of flight computer (RESTRT) 

(2) State estimation using sensor outputs (FILTER) 

(3) Guidance command calculation 

(a) Nominal plus pilot perturbation (PILOT) 

(b) Autopilot using prestored guidance 
trajectory (GUIDE) 

(4) Control command calculation 

(a) Proportional-Integral-Filter 
controller (CNTRL1) 

(b) Proportional-Integral 
controller (CNTRL2) 

(5) Control gain scheduling (GNADPT) 

(6) Filter gain scheduling (FLADPT) 

(7) State estimate extrapolation (FLPRP) 


1.7 TIME HISTORY SUMMARY 


Subroutine MSUMRY (Fig. 1.7-1) produces a printed 
state time history , (by calling HPRINT) and plots a time 
history of selected elements of continuous and discrete 
state variables (by calling . HPLOT) , 


R-31B16 



Figure 1.7-1 MSUMRY Block Diagram 
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PSUMRY (Fig. 1.7-2) prints the time history of selected 
elements of covariance matrices and plots selected covariance 
elements . 


ft— 31117 



Figure 1.7-2 PSUMRY Block Diagram 
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2 . 


TVHIS USAGE 


TVHIS may be run in one of two modes. A mean tra- 
jectory simulation examines the motions of the helicopter 
in a neighborhood of a reference flight path. Covariance 
calculation propagates the covariance from initial conditions 
along a specific flight path, assuming proportional-integral- 
filter control with velocity guidance. The program can be used 
to evaluate the design of filters, guidance algorithms, and 
digital-adaptive control laws. 


2.1 TVHIS INPUT 

Input to TVHIS is in two parts : card-image input on 

FORTRAN logical unit 5, and flight aerodynamic coefficients on 
FORTRAN logical unit 1. 

2.1.1 Card-Image Input 

Card-image data on logical unit 5 starts with a 
single NAMELIST ($STACK) read by the main program containing 
the number of stacked cases (NCASES). It is followed by 
NCASES blocks of data each containing 14 NAMELISTS and con- 
trol gain regression coefficients. 

There are only a few general rules which govern the 
specification of values for input parameters through a NAME- 
LIST. Letting STACK represent a general NAMELIST name, the 
first card in the input deck must contain a blank in the 
first column, $STACK starting in the second column, and a 
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blank in column 8. The specification of individual data 
items then begins in column 9 of the first data card and 
continues on succeeding cards if necessary. A data item 
for an unsubscripted variable is of the form v = c where 
v is the program mnemonic and c is a constant (e.g., NCASES=10 
and TIMRUN = 10.0). Data items for subscripted variables 
are of the form a = d^, ... , d^ or the form a(n) = d^ , ..., d m 

where a is the array name, n is an integer constant, and the 
d i are simple constants or repeated constants of the form k * c 
where k is the repetition , factor . With an unsubscripted array 
name, the total number of constants must equal the number of 
elements in the array (e.g., R = 5.0, 5.0, 15.0). With a 
subscripted array name, the total number of constants need 
not equal, but may not exceed, the number of array elements 
needed to fill the array starting with a(n) (e.g., TIMD(l) = 2.0 
and TIMG(8) = 1.E7, 2.5E5). All data items must be separated 
by commas, and any desired number of blanks may be inserted 
between data items. 

The end of a NAMELIST record is indicated by a $ after 
the last data item. The $ may appear directly after the last 
constant in the data item or may be separated by one or more 
blanks. For convenience in adding data items to a given NAME- 
LIST, a comma may always be inserted after the last data item 
and the $, signaling the end of the NAMELIST, may be placed on 
a single card. Furthermore, to make the end of the CDC NAME- 
LIST analogous to the end of the IBM NAMELIST, $END may be 
used in place of just a $. 

Three other important points regarding NAMELIST 
input data should be made. First, there are no restrictions 
on the order in which the data items are arranged within the 
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list. Second, the first column of every data card, not just 
the first card, must be left blank; all other 79 columns may 
be used. Third, the name of every NAMELIST to be read must 
appear in the input deck, even if no data items are specified 
for some NAMELIST. 

Assuming a multiple-case run, data input to override 
defaults for case N become defaults for case N+l and subsequent 
cases. For each of the NCASES runs, 14 namelists must 
appear in the input stream on logical unit 5, followed by 
control gain regression data. The format of this input deck 
is indicated in Fig. 2.1-1. NAMELIST input is described in 
Tables 2.1-1 through 2.1-21 specifying program mnemonic, a 
description of the variable, the units (if applicable) and a 
default value. 
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TABLE 2.1-1 


NUMBER OF STACKED CASES 
(NAMELIST/STACK/) 

T-1695 


PROGRAM MNEMONIC 

■ 

DESCRIPTION 

j DEFAULT 
™ ITS j VALUE 

NCASES 

| 

Number of stacked cases for this 1 
job 

- 

1 


TABLE 2.1-2 

INPUT /OUTPUT CONTROL 
(NAMELIST/NINOU/ ) 


T-1696 


PROGRAM MNEMONIC 

DESCRIPTION 

UNITS 

DEFAULT 

VALUE 

KIN 

Logical unit number for card-image 
input (if modified, subsequent 
card-image input will be read 
from logical unit KIN) 

- 

5 

KDAT 

Logical unit number for disk file 
containing aerocoef f icient data 

- 

1 

KOUT 

Logical unit number for printed 
output (if modified subsequent 
print statements will write to 
KOUT) 


6 

IGNPT 

Gain print flag: 

Produce reports of scheduled con- 
trol, filter gains if > 0 


1 

INOMPT 

Nominal state flag: 

Print time history of nominal 
state if > 0 

' 

1 

KPLOT[ 3] 

Plot control key 



KPLOT(l) 

Plots produced if > 0 

- 

0 

KPL0T(2 ) 

Not used 

- 

0 

KPLOT( 3 ) 

Not used 

- 

0 

KPRINT[2] 

Debug print control 



KPLOT(l) 

Not used 

- 

0 

KPLOT( 2 ) 

Produce filter debug printout if 
> 0 

— 

0 
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TABLE 2.1-3 

INTEGRATION CONTROL PARAMETERS 
(NAMELI ST/NRKINT / ) 


T-1697 


PROGRAM MNEMONIC 

DESCRIPTION 

UNITS 

DEFAULT 

VALUE 

ISTYPE 

Simulation type 

» 1 mean trajectory simulation 

Runge-Kutta-Gill integration 
“ 2 covariance calculation rec- 
tangular integration 
= 3 covariance calculation trape- 
zoidal integration 


1 

ACCU 

Accuracy limit for Runge-Kutta- 
Gill convergence test 

- 

0.0003 

DXIC[80] 

Initial value of state vector for 
Runge-Kutta-Gill integration 

- 

80*0.0 

WT [80] 

Vector of weights on state ele- 
ments for Runge-Kutta-Gill inte- 
gration 

NDIM 

Restriction: £ WT(I) = 1.0 

i=l 

where 


3*0.03, 

3*0.06, 

6*0.1621667, 

12*0.01, 

56*0.0 

i 

! 

• 

NDIM = 24 if actuator/rotor 
dynamics are on 
NDIM * 12 if actuator/rotor 
dynamics are off 
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TABLE 2.1-4 

WIND AND SENSOR FLAGS 
(NAMELIST /NREF / ) 

T-1698 


PROGRAM MNEMONIC 

| 

DESCRIPTION 

UNITS 

DEFAULT 

VALUE 

KTURB 

Wind/turbulence flag 
=0 no wind/turbulence 
= 1 wind/turbulence model 


0 

KSEN[3] 

Sensor key 

- 


KSEN(l) 

= 0 set observation vector 
ZVCOBS (31-46) to 0.0 
= 1 set ZVCOBS (31-46) to true 
state, control 


1 

KSEN( 2 ) 

Regulates angular rate , Euler 
angle sensors 
= 0 load true values 
= 1 corrupt true values with 
biases, noise 


0 

KSEN( 3 ) 

Regulates posit ion/ velocity 

sensors 

=0 no position/velocity sensors 

= 1 "perfect” accelerometer, 
altimeter, true airspeed, 
sideslip measurements 

= 2 corrupted accelerometer, air- 
speed, altimeter, sideslip 
measurements 

= 3 (1) above plus "perfect" MLS 

measurement 

* 4 (2) above plus noisy MLS 

measurement 

= 5 (not implemented) sensor 
dynamics 

= 6 (1) above plus "perfect" tri- 

lateration position, 
velocity 

= 7 (2) above plus noise trila- 

teration position, 
velocity 

i 

0 
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TABLE 2.1-5 

NOMINAL TRAJECTORY 
(NAMELIST/NMISS/ ) 


T-1699 


PROGRAM MNEMONIC 

DESCRIPTION 

UNITS 

DEFAULT 

VALUE 

XIND [3 , 42] 

Nominal position 


Table 

2.1-16 

XIND(1 , I) 

Axial position at waypoint I 

feet 

XIND( 2 , I) 

Lateral position at waypoint I 

feet 


XIND( 3 , I) 

Normal position at waypoint I 
(negative of altitude) 

feet 


VIND [3 , 42] 

Nominal velocity 


Table 

2.1-16 

VIND( 1 , I) 

Axial velocity at waypoint I 

feet /sec 

VIND( 2 f I) 

Lateral velocity at. waypoint I 

feet/sec 


VIND( 3 , I) 

Normal velocity at waypoint I 
(positive down) 

feet/sec 


TIMD[42] 

Scheduled waypoint times 


Table 

TIMD(I) 

Time at waypoint I 

secs 

2.1-16 

JTRAJ[42] 

Trajectory type 



• JTRAJ(I) 

Trajectory type following waypoint 




* 5 straight flight path 

- 



= 6 circular flight path 

- 



* 7 linear acceleration 
= 8 turn entry 

— 


NP 

Total number of waypoints 

- 

42 
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TABLE 2.1-6 

WIND INFORMATION 
(NAMELIST/NWIND/ ) 


T-1700 


PROGRAM MNEMONIC 

• DESCRIPTION 

UNITS 

DEFAULT 

VALUE 

TWIND[4] 

Array of time keys for wind inputs 


0 . 0 , 5.0, 

10.0, 15.0 

TWIND(I) 

Beginning time of I th wind input 
(relative to initial time of simu- 
lation ) 

secs 


TWIND(I + 1) 

Ending time of I th wind input 
(relative to initial time of 
simulation) 

secs 


VEWIND[3, 3] 

Array of wind inputs in earth axes 



VEWIHD(1, I) 

Axial component of wind between 
TWIND (I) and TWIND (I + 1) 

feet/sec 

0 0 
0 0 

H* 

O 

O 

VEWIND(2 , I) 

Lateral component of wind be- 
tween TWIND (I) and TWIND (I + 1) 

feet/sec 

i 

M 

O O 
O O 

O 

O 

VEWIND( 3 , I) 

Normal component of wind between 
TIMD (I) and TIMD (I + 1) 

feet/sec 

O 

O 

O 

O 

O 

QTURB [3,3] 

Turbulence driving noise covari- 
ance used for covariance calcula- 
tion 


9*0 . 0 


TABLE 2.1-7 




ROTOR /ACTUATOR DATA 
( NAMELI ST /NLIM / ) 


T-1701 

— 

PROGRAM MNEMONIC 

DESCRIPTION 

UNITS 

DEFAULT 

VALUE 

— 

LIMKEY[2] 

LIMKEY(l) 

Actuator limit flags 

0 *= off 

1 = on 

Actuator position limit key 


0 

— 

LIMKEY(2 ) 

Actuator rate limit key 

- 

0 


UMIN[4] 

Minimum actuator position vector 

inches 

- 2 . 0 , 0 . 0 , 


UMAX [4] 

Maximum actuator position vector 

inches 

-2.0, -2.0 
2.0, 4.0 


UDTLIM[4] 

UMAX, UMIN used when LIMKEY(1)«1 
Maximum absolute actuator rate 

inches/sec . 

2.0, 2.0 
3.0, 3.0, 


FROT [3] 

UDTLIM used when LIMKEY(2 ) = 1 
Rotor dynamics description array 


3.0, 3.0 
-576.0, 28. £ 


TACT [4] 

Actuator time constants 

secs 

1440.0 

4*0.0125 


KACT 

Actuator and rotor key 
=0 no rotor, actuator dynamics 
(state dimension NDIM « 12). 
= 1 rotor, actuator dynamics 
modeled (state dimension 
NDIM = 24) 


0 

i 
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TABLE 2.1-8 

INERTIAL CONSTANTS 
C NAMELIST /NINERT / ) 

(NOTE: ALL VARIABLES ARE REAL) 


T-1702 


PROGRAM MNEMONIC 

DESCRIPTION 

UNITS 

DEFAULT 

VALUE 

IX 

X-axis reference rotational iner- 
tia 

slug-f t 

37163.0 

IY 

Y-axis reference rotational iner- 
tia 

slug-ft 

20175.0 

12 

Z-axis reference rotational iner- 
tia 

slug-ft 

190059.0 

IXP | 

Present X-axis rotational inertia 

slug-ft 

37163.0 

IYP 

Present Y-axis rotational inertia 

slug-ft 

20175.0 

IZP 

Present Z-axis rotational inertia 

slug-ft 

190059.0 

IXZP 

Present X-Z product rotational 
inertia 

slug-ft^ 

201750.0 

MASS 

Reference mass 

slug 

1027.0 

MASSP 

Present mass 

slug 

1027.0 
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TABLE 2.1-9 

SENSOR DATA 
( NAMELIST/NSENS/ ) 


T-1703 


PROGRAM MNEMONIC 

DESCRIPTION 

UNITS 

DEFAULT 

VALUE 

XAACT[3] 

Location of MLS azimuth equipment 

MM 


XAACT(l) 

Axial position 


10000.0 

XAACT ( 2 ) 

Lateral position 


0.0 

XAACT ( 3 ) 

Normal position 

feet 

-12.0 

XEACT [3] 

Location of MLS elevation equip- 
ment 



XEACT(l) 

Axial position 

feet 

1000.0 

‘ XEACT(2) 

Lateral position 

feet 

300.0 

XEACT ( 3 ) 

Normal position 

feet 

-12.0 

XSACT 

Location of on-board accelero- 
meters forward of vehicle center 
of gravity 

feet 

0.83333 

ZSACT 

Location of on-board accelero- 
meters below vehicle center of 
gravity 

feet 

1.66667 

PHIB[3] 

- 

Diagonals of 3 * 3 accelerometer 
bias transition matrix 

■ 

3 * 9.99E-1 

PHIW [3] 

Diagonals of 3 x 3 wind transition 


3*9 . 96672E-1 
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TABLE 2.1-9 (CONTINUED) 

SENSOR DATA 
(NAMELIST/NSENS/) 


PROGRAM MNEMONIC 

DESCRIPTION 

UNITS 

DEFAULT 

VALUE 

q [22] 

Vector containing variance of 
non-time-varying sensor errors 
(biases , scale factors) 



Q(l> 

MLS range bias 

feet 2 

2.56E02 

Q(2) 

MLS azimuth bias 

o 

radians 

4 . 3681E-8 

Q(3) 

MLS elevation bias 

2 

radians 

1 . 2100E-4 

Q(4) 

MLS range scale factor 

- 

0.00E00 

Q(5) 

Accelerometer bias - body x-axis 

(feet/sec 2 ) 2 

6.40E-1 

Q(6) 

Accelerometer bias - body y-axis 

(feet/sec 2 ) 2 

6.40E-1 

Q(7) 

Accelerometer bias - body z-axis 

( feet/sec 2 ) 2 

6.40E-1 

Q(8) 

Accelerometer scale factor -body 
x-axis 

- 

6.15E-10 

QC9) 

Accelerometer scale factor - 
body y-axis 

- 

6.25E-10 

Q( 10 ) 

Accelerometer scale factor - 
body z-axis 

- 

6.25E-10 

Q(ll) 

Barometric altimeter bias 

feet 2 

1.0E-4 

Q(12 ) 

Barometric altimeter scale factor 

feet 2 | 

9 . OE-4 

Q(13) 

Body roll angular rate gyro bias 

(rad/sec) 2 

2 . 74576E-5 

Q( 14 ) 

Body pitch angular rate gyro bias 

(rad/sec) 2 

2 . 74576E-5 

Q( 15 ) 

Body yaw angular rate gyro bias 

(rad/sec) 2 

2 . 74576E-5 

Q( 16) 

Body roll angular rate gyro scale 
factor 

- 

4.00E-4 

0(17) 

Body pitch angular rate gyro 
scale factor 

- 

4.00E-4 

Q(18) 

Body yaw angular rate gyro scale 
factor 

- 

4.00E-4 

Q(19) 

Indicated air speed (IAS) bias 

( feet /sec ) 2 

2.50E-1 

Q(20 ) 

Indicated air speed (IAS) scale 
factor 

- 

2.50E-1 

0(21) 

Sideslip (3) bias 

radians 2 

2 . 74576E-5 

0(22) 

Sideslip scale factor 

- 

1.6E-2 
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TABLE 2.1-9 (CONTINUED) 

SENSOR DATA 
(NAMELIST/NSENS / ) 


T-1705 


PROGRAM MNEMONIC 

DESCRIPTION 

UNITS 

DEFAULT 

VALUE 

R[22] 

Vector containing variance of 
sensor noise* 



R(l)** 

MLS range or trilateration axial 
position 

feet 2 

0.0 

R(2) 

MLS azimuth or trilateration 

radians 2 or 

0.0 


lateral position 

feet 2 


R(3) 

MLS elevation or trilateration 

radians 2 or 

0.0 


normal position 

feet 2 


R(4) 

Trilateration axial velocity 

(feet/sec) 2 

0.0 

R(5) 

Trilateration lateral velocity 

( feet /sec) 2 

0.0 

R(6) 

Trilateration normal velocity 

(feet/sec) 2 

o.o 

R(7) 

Accelerometer body x-axis 

(feet/sec 2 ) 2 

4.00E00 

R(8) 

Accelerometer body y-axis 

( feet/sec 2 ) 2 

4.00E00 

R(9) 

Accelerometer body z-axis 

(feet/sec 2 ) 2 

4.00E00 

R( 10) 

Barometric altimeter 

feet 2 

6.25E02 

R(ll) 

Roll Euler angle 0 

radians 2 

8 . 4681E-6 

R( 12 ) 

Pitch Euler angle 6 

radians 2 

8 . 4681E-6 

R(13) 

Yaw Euler angle ip 

radians 2 

8 . 4681E-4 

R(14) 

Body roll angular rate p 

(rad/sec) 2 

8 . 4681E-G 

R( 15 ) 

Body pitch angular rate q 

(rad/sec) 2 

8 . 4681E-6 

R(16) 

Body yaw angular rate r 

(rad/sec) 2 

8 . 4681E-6 

R(17) 

Indicated airspeed IAS 

( feet/sec)“ 

2.8561' 

R(18) 

Sideslip 6 

radians 2 

7 . 62 129E-5. 

R( 19-22) 

Control position 

inches 2 

4*0.0 


* Interpretation of R for mean trajectory simulation. 

**R(l-6) depend on range and are calculated by subroutine SENINT. 
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TABLE 2.1-9 (CONTINUED) 

SENSOR DATA 
(NAMELIST / NSENS / ) 


T-1706 


PROGRAM MNEMONIC 

DESCRIPTION 

UNITS 

DEFAULT 

VALUE 

R [ 22 ] 

Vector containing variance of 
sensor noise* 



R(l) 

Axial position 

feet 2 

2.25 

R(2) 

Lateral position 

feet 2 

2.25 

R(3) 

Normal position 

feet 2 

2.25 

R(4) 

Body x-axis velocity 

(feet/sec) 2 

0.12 

R(5) 

Body y-axis velocity 

(feet/sec) 2 

0.64 

R(6) 

Body z-axis velocity 

(feet/sec) 2 . 

0.25 

R(7) 

Body roll rate 

(rad/sec) 2 

8 . 4681E-6 

R(8) 

Body pitch rate 

(rad/sec) 2 

8 . 4681E-6 

R(9) 

Body yaw rate 

(rad/sec ) 2 

8 . 4681E-6 

R(10) 

Roll Euler angle 

9 

radians 

8 . 4681E-6 

R( 11) 

Pitch Euler angle 

radians 2 

8 . 4681E-6 

R(X2) 

Yaw Euler angle 

r> 

radians 

3 . 0625E-4 


♦Interpretation of R for covariance run. 
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TABLE 2.1-10 


FLIGHT COMPUTER CONTROL SWITCHES 
(NAMELIST/NCFWCH) 


■ 

PROGRAM MNEMONIC 

DESCRIPTION 

UNITS 

DEFAULT 

VALUE 

. — 

* NRESTR 

signals flight computer reset when 
1 (set to 0 after reset - should 
be set to 1 in each case of multi- 

- 

1 



pie- 

case run) 



* — ' 

NAUTO 

Guidance indicator 

_ 

1 



*= 0 

automatic guidance 





* 1 

pilot input 




NCTYPE 

Control type selector 

_ 

1 



* 1 

proportional integral filter 
control law 





= 2 

proportional integral control 
law 




NGTYPE 

guidance type selector 

- 

1 



= 1 

velocity guidance 





= 2 

attitude guidance 




KTIl[4] 

Filter selection and mode setting 




KFIL(l) 

Angular filter control 

- 

1 * 



< 0 

load true Euler angles, body 






angular rates as estimated 




. 

= 1 

bvpass filter - copy sensed 
Euler angles, angular rates 
into estimates 





= 2 

use angular filter with full 
trigonometric functions 





= 3 

use angular filter with 
second-order approximations 
to sine, cosine functions 




KFIL( 2 ) 

Controls enroute filter 

- 

2 

- — ■ 


< 0 

bypass enroute filter 


1 



= 0 

’’perfect" enroute filter 
states-true values as esti- 


1 




mates 





= 1 

enroute filter, full trigo- 
nometric functions, sideslip 





- 2 

enroute filter, full trigo- 
nometric functions, zero 
sideslip 



! 


=3 

enroute filter, sideslip, 
second-order approximate 
trig functions 


i 



=4 

enroute filter , pseudo-zero 



— ■ 



sideslip, second-order ap- 
proximate trig function 




KFIL( 3 ) 

Controls trilateration filter 

_ 

2 



< 0 

bypass trilateration filter 





= 0 

true states as estimates 





« 1 

filtered states, full trigo- 
nometric functions, sideslip 





- 2 

filtered states, full trig 
functions, zero sideslip 





= 3 

filtered states, sideslip, 
second-order approximate 
trig functions 





* 4 

filtered states, zero side- 



— 

• 


slip, second-order approxi- 
mate trig functions 




KFIL(4 ) 

Controls MLS filter. Switch 
settings same as for KFIL(3) - 

- 

2 



trilateration filter 
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TABLE 2.1-11 

TIME PARAMETERS 
(NAMELIST/NTINT/ ) 


T-170H 


PROGRAM MNEMONIC 

DESCRIPTION 

UNITS 

DEFAULT 

VALUE 

DTGUID 

Time interval between successive 
guidance calls (PILOT or GUIDE) 

seconds 

0.1 

DTFIL 

Time interval between successive 
filter calls (measurement data 
rate) 

seconds 

0.1 

DTCNTR 

Time interval between successive 
calls to controller subrputine 

seconds 

0.1 

DTADPT 

Time interval between successive 
control gain calculations 

seconds 

1.0 

DFADPT 

Time interval between successive 
filter gain calculations 

seconds 

1.0 

DTOUTP 

Sample time interval for time 
history print /plot 

seconds 

0.05 

TIMRUN 

Elapsed time for simulation 

seconds 

1.0 

NWPF 

Waypoint number of nominal 
trajectory for start of simu- 
lation 


1 

TIMDEL ! 

Time delay between control com- 
mand calculation and its appli- 
cation 

seconds 

o 

o 

DTINT ' 

Starting integration step size 

seconds 

0.05 

TEVENT[6] 

Time key for pilot input 
EVMAG(I) will be applied for all 
TIME > TEVENT(I) 

seconds 

6*1 . 0E10 

ityp[s] 

Control command element which 
EVMAG is to perturb (1-4) 

- 

6*1 

EVMAG[6] 

j 

Magnitude of pilot event-will be 
applied as perturbation on nominal 
command 


6*0.0 
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TABLE 2.1-12 

PIF CONTROL COMMAND LIMITS 
(NAMELI ST/NCNTRL/ ) 


T-1709 


PROGRAM MNEMONIC 

DESCRIPTION 

UNITS 

DEFAULT 

VALUE 

COMLIM [2 , 4 ] 

Limits on control command output 
from PIF control law sent to 
actuators 



COMLI.M ( 1 , 1) 

Lower limit on first element of 
control UVCCOM(l) 

inches 

-5.0 

COMLIM ( 1 , 2) 

1 

1 

Lower limit on second element of 
control command UVCC0M(2) 

inches 

0.0 

COMLIM ( 1 , 3) 

Lower limit on UVCC0M(3) 

inches 

-5.0 

COMLIM ( 1 . 4) 

Lower limit on UVCC0M(4) 

inches 

-5.0 

COMLIM( 2 , 1) 

Upper limit on UVCCOM(l) 

inches 

5.0 

COMLIM( 2 , 2) 

Upper limit on UVCC0M(2) 

inches 

10.0 

C0MLIM(2 , 3) 

Upper limit on UVCC0M(3) 

inches 

5 . 0 

COMLIM ( 2 , 4) 

Upper limit on UVCC0M(4) 

inches 

5.0 
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TABLE 2.1-15 

COVARIANCE MATRIX INITIALIZATION 
(NAMELIST/NCOVAR/ ) 


T-1712 


PROGRAM MNEMONIC 

DESCRIPTION 

UNITS 

DEFAULT 

VALUE 

Pll[l2, 12] 

Helicopter rigid-body state 
covariance 

- 

PU(I,J)-0.0 
for lj«j 
P11(I,I)-1.0 

P12 [12 , ,5j 

Helicopter state-turbulence 
covariance 

- 

60*0.0 

P13 [12 , 4] 

Helicopter state-control covari- 
ance 

- 

48*0.0 

P14 [12 , 12] 

Helicopter state-measurement 
covariance 

- 

144*0.0 

P15 [12 , 4] 

Helicopter state-controller 
state covariance 

- 

48*0.0 

P22[5, 5] 

Turbulence covariance 

— 

P22(I , J)*0 . 0 
for I^J 
P22( I , I )-l . 0 

P23 [5 , 4] 

Turbulence-control covariance 

- 

20*0.0 

P24[5, 12] 

Turbulence-state measurement 
covariance 

- 

60*0.0 

P25[5, 4] 

Turbulence-controller state 
covariance 

- 

20*0.0 

P33[4, 4j 

Control covariance 

- 

16*0.0 

P34 [4 , 12] 

Control-state measurement 
covariance 

- 

48*0.0 

P35 [4 , 4] 

Control-controller state 
covariance 


16*0.0 

P44 [12 , 12] 

State measurement covariance 

- 

144*0.0 

P45 [l2 , 4] 

State measurement-controller 
covariance 

- 

48*0.0 

P55 f 4 , 4] 

Controller state covariance 

- 

16*0.0 
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TABLE 2.1-16 

HELICOPTER NOMINAL APPROACH TRAJECTORY (DEFAULT) 

T-1713 


WAYPOINT 

NUMBER 

TIME 

POSITION (FEET) 

VELOCITY (FEET/SEC) 

TRAJECTORY 

TYPE 

(i) 

TIMD(I) 

XIND( 1 , J) 

XIND(2,I) 

XIND( 3 , I ) 

VIND(1 , I ) 

VIND(2,I) 

VIND( 3,1) 

JTRAJ(I) 

1 

0.0 

-14000.00 

-400.00 

-750.00 

-77.78 

77.78 

0.0 

5 

2 

42.22 

-17283.16 

-716.50 

-750.00 

-77.78 

77.78 

0.0 

5 

3 

43.21 

-17360.94 

-638.72 

-750.00 

-77.78 

77.78 

0.0 

8 

4 

44.21 

-17439.39 

-561.61 

-750.00 

-79.80 

75.76 

0.0 

6 

5 

74.00 

-20400.00 

-600.00 

-750.00 

-77.78 

-77.78 

0.0 

6 

6 

104.30 

-20400.00 

-3600.00 

-750.00 

77.78 . 

-77.78 

0.0 

6 

7 

134.59 

-17400.00 

-3600.00 

-750.00 

77.78 

77.78 

0.0 

6 

8 

164.88 

-17400.00 

-600.00 

-750.00 

-77.78 

77.78 

0.0 

6 

9 

195.17 

-20400.00 

-600.00 

-750.00 

-77.78 

-77.78 

0.0 

6 

10 

225.47 

-20400.00 

-3600.00 

-750.00 

77.78 

-77.78 

0.0 

6 

11 

255.76 

-17400.00 

-3600.00 

-750.00 

77.78 

77.78 

0.0 

6 

12-34 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.0 

5 

35 

556.24 

-3206.47 

-0.7380 

-482.30 

110.00 

2.95 

8.333 

8 

36 

557.24 

-3096.47 

0.2460: 

-473.97 

110.00 

0.00 

8.333 

; . 5 

37 

557.62 

-3055.00 

0.00 

-470.83 

110.00 

0.00 

8.333 

7 

38 

558 , 62 

-2945.00 

0.00 

-462.84 

110.00 

0.00 

7.640 

7 

39 

592.67 

-403.05 

o.oo ! 

-202.80 

40.33 

0.00 

7.640 

7 

40. 

612.17 

-0.25 

0.00 

-53.82 

0.01 

0.00 

7.640 

7 

41 

613.17 

0.00 

0.00 

-50.00 

0.00 

0.00 

0.00 

5 

42 

700.00 

0.00 

0.00 J 

-50.00 

0.00 

0.00 

0.00 

5 


Note: Waypoints 12-34 are not specified: default is 0.0 
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TABLE 2.1-17 

GUIDANCE HELICOPTER APPROACH TRAJECTORY (DEFAULTS) 


WAYPOINT 

NUMBER 

TIME 

POSITION ( FEET ) 

VELOCITY ( FEET / SEC ) 

( I ) 

TIMG ( I ) 

XIG ( 1 , I ) 

XIG ( 2 , I ) 

XIG ( 3 , I ) 

VIG ( 1 , I ) 

VIG ( 2 , I ) 

VI G ( 3 , I ) 

1 

0.00 

- 14000.00 

- 4000.00 

- 750.00 

- 77.78 

77.78 

0.00 

2 

43.71 

- 14000.00 

- 600.00 

- 750.00 

- 77.78 

77.78 

0.00 

3 

74.00 

- 20400.00 

- 600.00 

- 750.00 

- 77.78 

77.78 

0.00 

4 

104.30 

- 20400.00 

- 3600.00 

- 750.00 

- 77.78 

- 77.78 

0.00 

5 

134.59 

- 17400.00 

- 3600.00 

- 750.00 

77.78 

77.78 

0.00 

6 

164.88 

- 17400.00 

- 600.00 

- 750.00 

- 77.78 

77.78 

0.00 

7 

195.17 

- 20400.00 

- 600.00 

- 750.00 

- 77.78 

- 77.78 

0.00 

8 

225.47 

- 20400.00 

- 3600.00 

- 750.00 

77.78 

- 77.78 

0.00 

9 

255.76 

- 17400.00 

- 3600.00 

- 750.00 

77.78 

77.78 

. 0.00 

10 

286.05 

- 17400.00 

- 600.00 

- 750.00 

- 77.78 

77.78 

0.00 

11 

306.62 

- 19000.00 

1000.00 

- 750.00 

- 77.78 

77.78 

0.00 

12 

310.89 

- 19292.89 

1292.89 

- 750.00 

- 77.78 

77.78 

0.00 

13 

327.62 

- 2000.00 

3000.00 

- 750.00 

0.00 

110.00 

0.00 

14 

370.46 

- 17000.00 

6000.00 

- 750.00 

110.00 

0.00 

0.00 

15 

379.55 

- 16000.00 

6000.00 

- 750.00 

110.00 

0.00 

0.00 

16 

388.65 

- 15000.00 

6000.00 

- 750.00 

110.00 

0.00 

0.00 

17 - 

417.21 

- 13000.00 

4000.00 

- 750.00 

0.00 

- 110.00 

0.00 

18 

435.39 

- 13000.00 

2000.00 

- 750.00 

0.00 

- 110.00 

0.00 

19 

452.62 

- 12292.89 

292.89 

- 750.00 

77.78 

- 77.78 

0.00 

20 

456.39 

- 12000.00 

0.00 

- 750.00 

. 77.78 

- 77.78 

0.00 

21 

482.10 

- 10000.00 

- 2000.00 

- 750.00 

77.78 

- 77.78 

0.00 

22 

485.87 

- 9707.11 

- 2292.89 

- 750.00 

77.78 

- 77.78 

0.00 

23 

503.11 

- 8000.00 

- 3000.00 

- 750.00 

110.00 

0.00 

0.00 

24 

520.34 

- 6292.887 

- 2292.89 

- 750.00 

77.78 

77.78 

0.00 

25 

524.11 

- 6000.00 

- 2000.00 

- 750.00 

77.78 

77.78 

8.333 

26 

542.11 

- 4600.00 

- 600.00 

- 600.00 

77.78 

77.78 

8.333 

27 

556.74 

- 3151.47 

0.00 

- 478.13 

110.00 

0.00 

7.640 

28 

558.12 

- 3000.00 

0.00 

- 50.00 

0.00 

0.00 

0.00 

29 

612.67 

0.00 

0.00 

- 50.00 

0.00 

0.00 

0.00 

30 

700.00 

0.00 

0.00 

- 50.00 

0.00 

0.00 

0.00 
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TABLE 2.1-18 
ANGULAR FILTER GAINS 

T-1715 


CONSTF (1-3) 

CONSTF (4-9) 

ANGULAR RATE GAINS 

EULER ANGLE GAINS 

0 . 96500E00 

0 . 27000E00 

0 . 97200E00 

0 . 2 7000E00 

0 . 95400E00 

0. 51000E-1 


0. 10000E00 


0 . 10000E00 


0. 10000E00 


TABLE 2.1-19 

ENROUTE FILTER GAIN REGRESSION COEFFICIENTS 


- T-1716 


INTERCEPT 
CONSTF (10-20) 

z 

RCER (1-11, 1) 

1(1 + C*VSQ) 
RCER (1-11, 2) 

9 . 31 774E-01 

0. 1 491 25-05 

— 0.96219 E— 0 2 

C . 49325E-02 

0. 38229E— 06 

-0. 18865E-02 

-0.31335E-0 3 

-9.391 26E-07 

C . 1 2 59 8 E— 03 

0. 9491 9E-01 

-0. 5C738E-0 8 

-0. 19297E-02 

0. 96315E-01 

-0. 35453E-07 

-0.17317E— 01 

Dt 1 0033F 90 

0.1 2026E-06 

0 .39365E-01 

6 . 623291-01 

0.72011 E— 07 

9.34871E 00 

-0 # 1 432 2 E — 01 

-0. 1 9 639 E- 07 

-0 .661 74E-02 

-0 . 7 3 72 6E — n 2 

-0.12211 E— 0 6 

-0 .59000E-01 

0. 4 995 0 E — 0 2 

0.0 

0. 0 

0. 999 ^ 9 E - 9 1 

0,0 

0.0 
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TABLE 2.1-20 

TR I LATE RAT I ON FILTER GAIN REGRESSION COEFFICIENTS 


INTERCEPT 
CONSTF (21-54) 

1/(1 + (C*R)**E ) 
RCTRI (1-34, 1) 

1/(1 + C*R ) 
RCTRI (1-34, 2) 

0. 87 26 6E — 0 2 

0 *551803 — 01 

i 

0 . 0 

0 . 20945E-02 

0.0 

0.33055E-01 

0 . 30S86E — 01 

-0 • 349 6 0 E — 0 2 

0 . 0 

0*9491 65-01 

0.0 

— 0 . 44 1 78E— 0 1 

O . 9461 3E-0 1 

n # r\ 

-0.41 878E-01 

0. 2S455F -02 

0. 87670E-01 

0. 0 

0. 0 

0.0 

0.0 

0. 0 

0.0 

0.0 

-0. 400C0E-03 

0.0 

0. 0 

-0. 35000E-03 

0.0 

0.0 

-0 . 2 860 OE— 02 

0.0 

0.0 

— 0 *2860 OE 00 

0.0 

0.0 

-0. 1 6 55 OF 00 

-0 .72509E-01 

0.0 

-0 . 30004E-C3 

0.0 

0.0 

0 • 73 0 0 OE 00 

0.0 

0.0 

0* 8300 OE 00 

0.0 

0.0 

0 . 3S425E — r 3 

0. 17883E-03 

0.0 

0 . 37 3 5 3E — 03 

0.0 

-0. 23919E-03 

0. 4 3 86 OE — 0 2 

0.0 

0 . 2344 4 E— 0 2 


0.0 

0 . 82652 E 00 

C. 76743E-0 1 

0.0 

0.82647E 00 

D. 1 0 0 0 OE — 0 3 

0 • 13964E on 

0.0 

0 • 39584E-02 

0. 0 

-0. 39584E-02 

0 . 46956E-02 

0.0 

-0.46956E-02 

-0 . 1 6274E-04 

-0 . 45481 E-04 

0.0 

-0. 1321 2E-04 

0. 0 

0.0 

-0 • 2320 7E-03 

0.0 

-C . 41 808E-03 

-0 . 1 O239E-01 

0.0 

-0.14056E 00 

-0. 10229E-01 

0 • 0 

-C. 14056E 00 

-0* 1 00 COE— 04 

-0.1 949CE-0! 

0 .0 

-0. 66239E-03 

0.0 

0 . 662 39E— 0 3 

-0. 7921 6E-03 

0.0 

0. 7921 6E-02 

0 . 4 9 95 0 E— 0 2 

0 

0 .0 

0. 99900E-0 1 

0.0 

0.0 
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TABLE 2.1-21 

MLS FILTER GAIN REGRESSION COEFFICIENTS 


INTERCEPT 
CONSTF (55-109) 

1/(1 + (C/10000 )**2 ) 
RCMLS (1-55, 1) 

1/(1 + (C/4000 )**2 ) 
RCMLS (1-55, 2) 

0 • 29967E— 0 1 

-0 .299675-01 

0.0 

0 . 50 43 3 5 — 0 1 

0.0 

0 . 241875-02 

0. 0 

0.0 

0.0 

o.o 

0.0 

0.0 

n • O 

0.0 

0.0 

0. 841 1 05-0 1 

0.0 

0.92004E-01 

0.0 

0.0 

0.0 

0 . 0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.547535-01 

0.0 

0.40720E 00 

0. 554765-01 

0.0 

-0. 523875-01 

0. 0 

0.0 

0.0 

O.o 

0.0 

0.0 

0.211835 00 

— 0 • 2 0393E 00 

0.0 

-0. 1 2 83 0 E — 0 2 

0. 1 2 839 E— 02 

0.0 

-0. 13687E-01 

n • 0 

0 . 77826 E— 02 

0. 0 

0.0 

0.0 

0.0 

0.0 

0.0 

0 . 0 

0 . 0 

0.0 

0. 1451 55 00 

— 0 • 20 309 E 00 

0.0 

0.0 

0.0 

o.o 

0.0 

0.0 

0.0 

0. 0 

0.0 

0.0 

-0.1 1 4 0 0 E — 0 2 

0.0 

-0 . 13700 E— C 1 

0.262395 00 

0.0 

-0.24156E 00 

0 . 0 

' 0.0 

0 . 0 

0 . 0 

0.0 

0.0 

0.402075 00 

0.0 

— 0 . 320 34 E 00 

0. 625025-02 

-0 . 62502E-02 

0.0 

0.131 755-01 

0.0 

0. 131 08E-02 

0. 0 

0.0 

0. 0 

0 . 0 

0.0 

0.0 

0.0 

0.0 

— 0.1 4500 E 01 

0. 3741 25-0 1 

0. 0 

0. 1 1 095E 00 

O.o 

0.0 

-0.1 59 3 0 E— 0 2 

-0.50568E 01 

0 . 0 

-0.94016E 02 

0.0 

0.0 

-0. 16427E-02 

0.0 

0.0 

0 .817655 00 

0. 578575-01 

0.0 

-0. 55 1 74 E— 0 1 

0. 0 

0.0 

0.0 

0 . 0 

0.0 

0.0 

0.661585-01 

0.0 

-0. 61397E-01 

-0. 4Q51 6 E — 0 3 

0 .48516E-03 

0.0 

-0.117! 2E-02 

0.0 

-0 . 152745-03 

0. 0 

0.0 

0.0 

0 . 0 

0.0 

0 .0 

0.11 3845-02 

0.0 

0.27304E 00 

-0. 4276 15-02 

0.0 

-0. 15984E-01 

0 . 0 

0.0 

0 . 24 86 6 E— 0 3 

-0 . 341485-01 

0.0 

0.105875 02 

0. 0 

0. 0 

0. 25607E-03 

0 . 0 

0.0 

-0.125985 00 

-0. 820475-02 

0.0 

0.783795-02 

0. 0 

0.0 

C .0 

0.0 

0.0 

o.o 

-0 . 1 0 2 3 3 E— 0 1 

0.0 

0 . 93 1 69 E— 0 2 

0. 49950E-02 

0.0 

o.o 

0 . 999005-0 l 

0.0 

0.0 
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Following all NAMELIST input, subroutine RDDAT 
reads values for control gain regressions into arrays 
REGCOF and CONST. Values for REGCOF are read 
row-wise, followed by values for CONST. If NCTYPE = 1 
(PIF controller) 84 rows for REGCOF are read, followed 
by 84 values for CONSTF. If NCTYPE = 2 (PI controller) 

52 rows for REGCOF are read', followed by 52 values for 
CONSTF. Each row of REGCOF is input using 3 card-image 
records, as follows: 

Record 1 columns 1-10 must be blank or numeric - 
(110). -The information contained 
therein is ignored. 

Record 2 REGC0F(I,J) J = 1, 2, 3, 4 are read 
under the FORMAT: E12.6, 6X, E12.6, 

6X , E12.6, 6X , E12.6, 6X. 

Record 3 REGCOF ( I, J) J = 5, 6, 7, 8 are read 
under the FORMAT: E12.6, 6X, E12.6, 

6X , E12.6, 6X , E12.6, 6X. 

After 84 or 52 rows of REGCOF have been read (252 
or 156 card-image records), values for CONSTF(I) are input, 

4 to a record, under FORMAT: E12.6, 6X, E12.6, 6X, E12.6, 

6X. CONSTF(I) contains the regression constant for the 
I^* 1 gain. REGCOF(I,J) contains the sensitivity of the I*"* 1 
gain the the scheduling variables. The scheduling 

variables are: 

1. VT-Total air-relative velocity (ft/sec) 

2. W-Body z-axis velocity (ft /sec) 

3. vt** 2 Square of air-relative velocity ( (ft/sec)**2) 

4 . PSIDOT-Body yaw Euler angle rate of change 
(rad/sec) 

5. VT**2 * PSIDOT 
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6. W * PSIDOT 

7. VT**2/(1 + C*VT**2) 

8. 1/(1 + C*VT**2) 


where : 

C = 0.123E-3 

2.1.2 Permanent Disk File 

The disk file referenced by logical unit 1 contains 
data in a core-image format. It is read by subroutine FLTINT 
using an unformattted FORTRAN READ statement. The file con- 
tains 297 different cases of aerodynamic coefficients. Each 
case consists of 71 aerodynamic coefficients used to deter- 
mine stability derivatives for the system matrix F. Each of 
the 297 different cases corresponds to a specific flight 
condition. The 3 variables which determine the exact flight 
condition are: 

1. i't -Rate of heading change (rad/sec) 

2. VH-Horizontal velocity (knots) 

3. VV-Vertical velocity (feet /min) 

FLTINT interpolates between the aerodynamic coefficients 

for the two flight conditions which bracket the current flight 

condition. 


2.2 TVHIS OUTPUT 

Output from TVHIS is written to FORTRAN logical 
unit 6, unless variable KOUT has been changed in the namelist 
NINOU in which case all output is written to logical unit KOUT.. 
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Output is functionally divided into 5 types: 

1. Input data base summary from RDDAT 

2. Special event reports from specific subroutines 

3. Time history print 

4. Time history plot 

5. Error messages. 

Each of these output types is described below. 

2.2.1 Input Data Base Summary 

Subroutine RDDAT produces a report on the TVHIS data 
base at the beginning of each case. Description paragraph 
numbers correspond to numbered areas in the accompanying figures 

Simulation control variables are printed as shown in 
Fig. 2.2-1. 

1. The first data printed delimits the scope of the 
flight segment being considered; the initial waypoint on the 
nominal trajectory, and the final time along the nominal tra- 
jectory are printed. 

2. The integration state dimension, NDIM, is printed 
next. For covariance calculation (ISTYPE = 2 or ISTYPE = 3) 

NDIM = 500. For a mean trajectory simulation, (ISTYPE = 1) 

NDIM depends on flag KACT from NAMELIST/NLIM/ . For KACT = 0, 
NDIM = 12 (12 rigid-body states will be used). For KACT = 1, 
NDIM =24 (12 rigid-body states, 4 actuator states, 8 rotor 
states). 
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INPUT RE PD PT FROM RDDAT 


1. FLIGHT FROM WAYPOINT NO. 1 TO A TRAJECTORY TIME OF 2.0 SECS 

2. DIMENSION OF STATE: 500 

3a. COVARIANCE CALCULATION: ISTYPE = 2 

CONTINUOUS STATES 

12 RIGID BODY STATES 
5 TURBULENCE STATES 

DISCRETE STATES 

4 CONTRPL COMMANDS 

to ' 12 PIGID-BODY STATE ESTIMATES 

<0 4 INTEGRATOR STATES 

3b. MEAN TRAJECTORY SIMULATION: ISTYPE = 1 

INTEGRATION ACCURACY LIMIT IS 0 .30000 POE-03 
INTEGRATION STATE WEIGHTS 

0. 3235246E-02 0 . 3 236246E- 02 0 . 3236246E-02 0 . 6472491E-02 0 . 647 249 1 E- 0 2 

^.64724915-02 0.1618212E 00 0.1618212E 00 0.1618212E 00 0.1618212E 00 

0.1618212E 00 0.1618212E 00 

4, OUTPUT TIME INTERVAL - DTOUTP = 0.05000 SECS 

Figure 2.2-1 Simulation Control Variables 
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3a. For a covariance calculation, ISTYPE indicates 
the integration method. 

ISTYPE = 2 Rectangular Integration 

ISTYPE = 3 Trapezoidal Integration 

3b. For a mean trajectory simulation, the accuracy 
limit, ACCU for Runge-Kutta integration, as well as state 
weighting vector WT, are printed. 

4. Output time interval, DTOUPT, is the time between 
successive samples for time history PRINT/PLOT. 

Pilot input time history and the nominal trajectory 
are printed in the format of Fig. 2.2-2. 

5. A pilot perturbation command history is printed 
if guidance is manual and ISTYPE = 1 

COMMAND VECTOR ELEMENT = ITYPE(I) 

TIME = TEVENT(I) 

PERTURBATION MAGNITUDE = EVMAG(I) 

The perturbation magnitude EVMAG(I) will be applied 
as a perturbation on the ITYPCI)^* 1 element of the nominal 
guidance command for all TIME >, TEVENT(I). 
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rMLW PFP T UP HAT I UN COMMAND HI5TOPY 


COMMAND 

VECTOR 

element = 3 

TIME = 

o.iooooe ii 

PERTUnOAT 1 ON 

MAGNI TURE 

COMMAND 

VECTOR 

element = I 

TIME = 

0 . 1 0000E 1 1 

PERT URD AT 1 ON 

MAGHI T UDE 

COMMAND 

VECTOR 

ELEMENT =- 1 

T | ME - 

O.IOOOOE 11 

PFpTUnnATl ON 

MAGNI 1 ODE 

COMMAND 

VECTOR 

ELEMENT - 1 

TIME = 

0. 1 OOOOE 1 1 

PER TUROATI ON 

MAGNI TUDE 

COMMAND 

VECTOR 

ELEMENT - 1 

TIME = 

0. 1 OOOOE 1 1 

PER TUROAT ION 

MAGNI TUDE 

COMMAND 

VECTOR 

element = 1 

TIME = 

O. 1 00 n 0E 1 1 

PERTUROAT ION 

MAGNI TUDE 


0.0 
0.0 
0 .0 
0.0 
o .n 
0 . 0 


6 . NOMINAL TPAJFCTOPY 

tO WAYPOINT NO. I 
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T I Mr - 

0. 0 

SEC 5 

POSITION » 
VELOCITY = 
TRAJECTORY 

-0* I A OOOOOE 05 
-0.77T0PO0E n 2 
type: straight 

-O-AOOOOOOE 0 A 
0.7776O00E C 2 
tlight JTRAJ - 

-0. 7500000E 
0 . 0 
5 

03 

EEET 
r EET 

PEP 

SEC 


WAYPQ I Ml 
TIM'’: r 

NO. 2 

* 2.21 

SECS 

POSITION = 
VELOCITY = 
trajectory 

-0 . 1 720316E 05 
-0 .77 7800QE 02 

type: straight 

-0.71 G5000E 03 
0.777 00 OOE 02 
TLIGHT JTPAJ - 

-0. 750 OOOOE 
0.0 

5 

03 

FEET 
EEC T 

PER 

SEC 


W AYPD I tJI 
t 1 M r - 

NO. 3 
A3. 21 

r.^cs 

P0S1 T ION = 
VELOCITY = 

-0.173609AE 05 
-0.7770000E 02 

— 0.6307200 E 03 
0. 77 70000E 02 

-0. 75 0000 OE 
0.0 

03 

FEET 

FEET 

PEP 

SEC 


tpajfctopy iypf: turn entry jtraj - n 


Figure 2.2-2 


Pilot Input 


Nominal Trajectory 
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6. NOMINAL TRAJECTORY is printed in all cases, 

(Fig. 2.2-2) starting with the initial waypoint, and con- 
tinuing with all waypoints in the nominal trajectory. 

WAYPOINT NUMBER I is printed, followed by: 

TIME = TIMD(I) 

EARTH-RELATIVE POSITION = XIND(1,I), XIND(2,I), XIND(3,I) 

EARTH-RELATIVE VELOCITY = VIND(1,I), VIND(2,I), VIND(3,I) 

TRAJECTORY TIME = trajectory type of the 

segment between waypoint I 
and waypoint 1+1. 

If the current case is a covariance calculation, then 
noise matrices are printed, as described in Fig. 2.2-3. 

7. Diagonals of turbulence driving noise, QTURB(3,3) 
are printed. 

8. Variance of sensor noise on 12 rigid-body states 
R(l-12) is printed. Position noise R(l-3) affects earth- 
relative position x, y, z. Velocity noise R(4-6) affects 
body-axis velocity u, v, w. Angular rate noise R(7-9) af- 
fects body angular rates p, q, r. 

Angular noise R(10-12) affects Euler angles 4> , 0, ip . 

RDDAT will then call COVPT to print the normalized forms of 
initial covariance matrices. 

If the current case is a mean trajectory simulation 
aircraft and sensor data are printed in the format of Fig. 2.2-4. 

9. If no wind inputs are specified (KTURB = 0) this 
fact is indicated. If KTURB = 1, scheduled wind inputs are 
printed. 
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7. 


TURBULENCE DRIVING NCISE VARIANCE 


0 T U R B ( 1.1) = 

OTUR B ( 2 • 2 ) = 

QTRU B (3*3) = 


0 . 0 
9.0 
0.0 


8. EFFECTIVE NOISE ON FILTER OUTPUT : 

POSITION NOISE VARIANCE 

R ( 1 ) = 0.0 FT**2 

R ( 2 ) = 0.0 FT**2 

R ( 3 ) = 0.0 FT**2 

VELOCITY NOISE VARIANCE 

R ( 4 ) = 0.0 FPS**2 

P C 5 ) = 0.0 FPS**2 

R 1 6 ) = 0.0 FPS**2 

ANGULAR RATE NOISE VARIANCE 

R ( 7 ) = 0 • 400 0 00E 01 RAD/SEC**2 

RC8) = 0 . 4 000 00 E 01 RAD/SEC**2 

R ( 9 ) = 0.40 00 00 E 01 RAD/SEC**2 

ANGULAR NOISE VARIANCE 

R ( 1 0 ) = 0. 625000E 03 RAD**2 

P(ll) = 0.846810E-05 RAD**2 

R ( 1 2 ) = 0.846819E-05 RAD**2 


Figure 2.2-3 


Noise Covariance Matrices 
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9. 


NO WIND OR TURBULENCE INPUTS 


10 , 

11 . 


to 


AIRCRAFT DATA 

PERFECT ANGLE AND ANGULAR RATE MEASUREMENTS 
PERFECT POSITION AND VELOCITY MEASUREMENTS 


13. 

14. 


NO ACTUATOR POSITION LIMITS 


NO ACTUATOR RATE LIMITS 


Figure 2.2-4 Aircraft and Sensor Data 


I 


J ) J J .1 1.1} 1 1 ) | 
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10. If noisy angular rate and angular sensors are in 
use, the variance of sensor noise on these states is printed. 

11. If accelerometers are modeled, their location 
relative to vehicle center of gravity is printed. 

12. If MLS or trilateration sensors are in use, 
the location of MLS azimuth and elevation towers is printed. 

13. If actuator displacement limits are to be used, 
(LIMKEY(l) =1) they are printed; if they are not in use, 
this fact is indicated. 

14. If actuator rate limits are to be used, (LIMKEY(2) 
1) the limits are printed; if rate limits are not used, this 
fact is indicated. 


Whether a mean trajectory simulation or a covariance 
calculation is being processed, the status of flight computer 
control memory is printed in the format of Fig. 2.2-5 (para- 
graphs 15 through 20). 


15. Flight computer control switches 

are printed. 

NRESTR 

should be 1 to indicate flight computer 
reset at beginning of case 

NAUTO 

0 indicates automatic guidance 

1 indicates manual guidance (nominal plus 
pilot perturbation) 

NCTYPE 

CONTROL TYPE 
= 1 PIF controller 
= 2 PI controller 

• 

NGTYPE 

= 1 velocity guidance 
= 2 attitude guidance 


KFIL(l) 

controls angular filter 
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FLIGHT COMPUTER DATA 


15 


FLIGHT COMPUTER CONTROL SWITCHES 


NRESTR = 
NAUTO = 
NCTYPE = 
NGTYPE = 


1 

0 

1 

1 


KFIL(l) - 1. 

KF I L ( 2 ) = -1 

KF I L < 3 ) = -1 

KFILH) = 1 


to 

i 

03 

CD 


16. 


17. 


GUIDANCE TIME INTERVAL 

DTGUID = 

0. 10000 

SECS 


CONTROL TIME INTERVAL 

DTCNTR = 

0 . 10000 

SECS 


FILTER TIME INTERVAL 

- 

DTFIL = 

0 . 10000 

SECS 


GAIN ADAPTATION TIME 

INTERVAL - 

DTADPT = 

1 .00000 

SECS 


FILTER ADAPTATION INTERVAL 

DFADPT = 

1 . 00000 

SECS 


CONTROL APPLICATION 

TIME DELAY 

-TIMDEL = 

0.0 

SECS 


CONTROL COMMAND LOWER 

LIMITS : 

-5.0000 

1 

o 

• 

o 

5.0000 

-5.0000 

CONTROL COMMAND UPPER 

LIMITS : 

5.0000 

10.0000 

5. 0000 

5.0000 


Figure 2.2-5 Flight Computer Control Settings 
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KFIL(2) controls enroute filter 
KFIL(3) controls MLS filter 
KFIL(4) controls trilateration filter 

16. Flight computer time intervals are printed. 

DTGUID time between successive guidance calculations 

DTCNTR time between successive control calculations 

DTFIL time between filter calculations-data 

sample rate 

DTADPT time between successive control gain 
calculations 

DFADPT time between successive filter gain 
calculations 

TIMDEL delay between time control command is 

calculated and time it is sent to actuators 

i 

17. Lower and Upper limits on PIF control command 
push-down are printed. 

18. Control gain regression coefficients. 

CONSTF and REGCOF are printed (Fig. 2.2-5). CONST appears 
under the heading "Intercept", and REGCOF is printed column- 
wise, with the gain-scheduling variable for each column as 

a header. 

19. If a mean trajectory simulation is in progress, 
the filter gain regression coefficients for each of the filters 
in use (angular, enroute, MLS, trilateration) are printed 

in the format of Fig. 2.2-6. For the angular filter, all 
gains are constant. For enroute, MLS, and trilateration filters, 
regression constants from CONSTF are printed under the heading 
"Intercept", while regression coefficients RCER, RCMLS, and 
RCTRI are. printed in two columns. The column heading indi- 
cates which scheduling variable is involved. 
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CONTROL GAIN REGRESSION COEFFICIENTS 


to 

I 

CO 

00 


INTERCEPT 

CAPV 

W 

CAPV-SQ 

PSIDOT 

- 0 . 30848 E -01 

- 0 . 20237 E -04 

0 . 37452 E -04 

0 . 0 

0 .0 
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0.0 

0. 0 

0.0 
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0.0 • 

0 .0 

0 . 30907 E-F 1 

0 . 21 6 C 6 E -05 

- 0 . 49992 E -04 

0 . o 

0.0 

- 0 . 42676 E -02 

0 . 34963 E -04 

0 . 1 75 5 2 E - 04 

0 . 0 

0.0 

0 . 28526 E -02 

- 0. 17 4 09 E— 0 4 

“ 0.504 77 E “ 0 4 

0.0 

0.0 

- 0 . 76 3 82 E — ° 1 

0. 3 3 0 52 E— 0 3 

0.0 

“ 0 . 61299 E -06 

0 . 0 

- 0 . 13502 E-C 2 

0.2 871 6 E - 0 4 

0 . 48202 E -04 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

O.o 

0 . 54155 E 00 

- 0 . 10457 E -02 

- 0 . 241 57 E -02 

0. 0 

0.0 

0.0 

0.0 

0.0 

0.0 

0 .0 

0.4358 7 E 00 

- 0 . 1 4592 E -0 1 

0 . 97721 E -03 

0.0 

0.0 

0 • 37 9 82 E 01 

0 . 1 0090 E -02 

0 . 22799 E -02 

0 . 96354 E -05 

“ 0 . 58140 E 

0 . 1 6926 E -02 

- 0 . 1 5 65 4 E - 05 

- 0 . 1 1 939 E -05 

0.0 

0.0 

0 . 180 1 5 E -0 2 

- 0 • ! 3856 E -04 

0 

0 . 24608 E -07 

0.0 


P S I DOT * VSQ 

PSI D0T*W 

VSQ/ ( 1 +C*\/SQ) 1/(1 +C*VSQ) 

-0 .79836E-06 

0.0 

-0. 98840E-06 

0 .0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0. 37048E-05 

0. 0 

-0 .50554E-06 

0.0 

0. 0 

-0. 19132E-02 

0.0 

“0.52693E-01 

0." 

0.0 

o.o 

0.0 

0. 1 0280E-05 

0.0 

-0.56398E-06 

0.0 

0. 2 1 26 2 E — 06 

0.0 

-0.1 2279E-05 

0.0 

0 . 1 81 30E-05 

0.0 

0 .64943E-06 

n.O 

0.0 

0 . 11 40 8E“ 0 3 

0.0 

0.30802E-01 

0.11 190E-04 

0.0 

-0 .641 32E-06 

0.0 

^.0 

0.0 

0.0 . 

0 .0 

0.3951 8E-03 

0.0 

“ 0 • 37 88 3E— 0 4 

0.0 

0.0 

0.0 

0.0 

0.0 

16762E-03 

0.0 

0 . 36 8 3 5E“ 0 3 

0.0 

0.0 

0.0 

0.0 

0.0 

0.11 553E-07 

0.0 

0 • 78379E-08 

0.0 

0 . 0 

-0.55982E-05 

0.0 

“0.17413E-02 


Figure 2.2-6 Control Gain Regression Coefficients 
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ENPOUTE FILTER 
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-0.207007 00 

0.0 


Figure 2.2-7 
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20. If guidance is automatic (NAUTO = 0), the full 
guidance trajectory is printed as in Fig. 2.2-7. Earth-rela- 
tive position and velocity vectors for the waypoint are 
followed by the 4x4 guidance gain matrix printed columnwise, 

2 columns per line. 

WAYPOINT NO = guidance waypoint I 

TIME = TIMG(I) 

POSITION = XIG( 1 , I ) , XIG( 2 , I ) , XIG(3,I) 

VELOCITY = VIG( 1 , I ) , VIG(2,I), VIG(3,I) 

2.2.2 Special Event Reports 

Special event reports from subroutines are described 
in Tables 2.2-1 through 2.2-11. The contents are: report 

title, description of variables in the report, and the sub- 
routine originating the report. 

2.2.3 Time History* Printout 

Time history printout appears as a series of pages, 
each with up to 9 columns of data. The first column is the 
sample time, and each succeeding column represents the value 
of a specific variable at the time given in column 1 (Fig. 2.2-8). 

2.2.4 Time History Plot 

Time history plots are printer plots with the time 
axis described down the page, and the values of up to 3 variables 
represented by separate symbols across the page. The distance 
between the maximum and minimum value of the variable is repre- 
sented by 120 print positions across the page, Fig. 2.2-9. 

Labels, descriptions, and program source variables for time 
history printout and plots are described in Table 2.2-12. 

Time history prints and plots occur either at the end of a 
run, or whenever 201 samples of program data have been taken 
and the intermediate storage arrays TIM, PXC, and PXD are full. 
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OEL-O 

OEL-C 

XE 

Z E 

THETA 

U 

Q 

W 

SEC 

IN 

( N 

f t 

FT 

OEG 

FPS 

DEG/SEC 

FPS 

0.0 

-0.7231 

3.2171 

-14000. 0000 

-750. 0000 

2.9380 

109. 8529 

0.0 
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3.21 71 
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0.0 
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0 . 1 000 
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Figure 2.2-8 
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TIME (SECS) O “ DETA-T ( OEG ) * = OETA-EST ( REG I 

-0.0584 -0.0477 -0.0370 -0.0263 -0.0156 -0.0049 0.0050 0.0165 0.0272 0.0379 0.0406 

0.0 f 4- 4- 4- — * 40 4 4 4 4 4 — f 


0.5000 


to 

to 1.0000 


I .5000 


2. 0000 


Figure 2.2-9 Sample Time History Plot 
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2.2.5 Error Messages 

There are three possible error messages from TVHIS. 

All indicate fatal errors, and all will be followed by abnormal 
program termination. The first error message is caused by an 
attempt to run a mean trajectory simulation with an invalid 
control-guidance combination. The three correct combinations 
of control type and guidance type are printed along with the 
entered control/guidance type, and program execution is termi- 
nated by RDDAT. 

A second error message is a result of an attempt to 
run a covariance calculation with NCTYPE f 1 or NGTYPE f 1. 

The covariance routines are only programmed for PIF velocity 
control. Program is terminated by RDDAT. 

A third error message is caused by failure of the 
Runge-Kutta-Gill integration routine RKGS to converge. Subrou- 
tine RKGS first integrates a state vector over a given interval 
using a step-size of DTINT and then repeats the integration 
using a step-size of DTINT/2. Letting SI and S2 represent the 
state vectors after the two integrations, respectively, the 
following sum is computed: 

NDIM 

E = l WT ( I ) (S2 ( I ) - S1(I) ) 

1 = 1 

If E > ACCU, the specified accuracy limit, the value of DTINT 
is halved, and the whole process is repeated. If this conver- 
gence test fails after 10 bisections of the step-size DTINT, 
RKGS will return a value of 11 in variable IHLF. If the 
main program receives a value of 11 for IHLF from RKGS, it 
will call MSUMRY for the final PRINT/PLOT for this case and 
then abnormally terminate the program, indicating that there 
was an error in RKGS . 
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Corrections for the first 2 types of errors involve 
entering the proper values of NCTYPE and NGTYPE in NAMELIST 
/NCFWCH/ . The third type of error may be due to one of the 
following: 

• Accuracy limit ACCU in NAMELIST/NRKGS/ 
was smaller than machine precision 
allows 

• State weight vector WT contained improper 
values (recall ][ WT(I) = 1.0) 

I 

• There was an abrupt discontinuity in the time 
history of one or more elements of the state 
vector possibly due to: 

- unreasonable perturbation EVMAG in NAMELIST 
/NTINT/ 

- incorrect control gains input 

- incorrect guidance gains GVD in NAMELIST 
/NGUIDE/ 

- incorrect filter gains in NAMELIST/NFILDT/ 

- unreasonable wind magnitude VEWIND in 
NAMELIST/NWIND/ . 
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TABLE 2 . 2-1 

SPECIAL EVENT REPORT - MI SI ON 


T-1694 


REPORT TITLE 

DESCRIPTION - VARIABLES 

SUBROUTINE 

DESCRIPTION OF 
WAYPOINT AT BE- 
GINNING/END OF 
SEGMENT 

TIME 

HORIZONTAL VELOCITY - VH - magnitude of velo- 
city in horizontal plane 

VERTICAL VELOCITY - W - velocity in vertical 

plane 

! 

HEADING - PSI - course or heading j 

of helicopter j 

HEADING RATE - PSIDOT - time rate of 

heading change 

HEADING ACCELERATION- PSIDOT - time rate of 

change heading rate j 

SYSTEM MATRIX F 12x12 system dynamics 1 

matrix 

INPUT MATRIX G 6x4 matrix determining ^ 

effect of control input ; 

on derivative of "inner 6" 
states 

TURBULENCE INPUT 6x3 matrix determining 

MATRIX HT effect of wind turbulence 

on derivatives of . 

"inner 6" states 

NOMINAL STATE - XVNO , XVN1 - nominal helicopter, 

state vector 

XVNO (1-3) Helicopter position - x, y, z 

XVNO (4-6) Body-axis velocity - u, v, w 

XVNO (7-9) Body angular rates - p, q, r 

XVNO (10-12) Euler Angles - <i> , 6, ty 
XVNO (13-16) Control Position - 6g, 6^ f 6g, 

XVNO (17-20) Control rate - 5g, 5 C , <5g, <5 R 

XVNO (21-24) Actuator command - 5g A , o^ A , ^SA’ 

£ ra 

DYNAMIC TRIM ANGULAR RATES -XDYNO (4-6) or 

XDYN1 (4-6) 

DYNAMIC TRIM EULER ANGLES -DNUSTO or DNUST1 

DYNAMIC TRIM CONTROL VECTOR - DUST0 or DUST1 

MIS ION 
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TABLE 2.2-5 

SPECIAL EVENT REPORT - FLADPT 

T-1723 


REPORT TITLE 

DESCRIPTION - VARIABLES 

SUBROUTINE 

REPORT FROM 

SUBROUTINE 

FLADPT 

ANGULAR FILTER GAINS FILGA 

FILGA(I) = CONSTF(I), I = 1, 9 
MLS FILTER GAINS FLGMLS 

FLGMLS ( I ) = CONSTF( 54 + I) 

+ (1 + (RE/100000) 2 )*RCMLS(1, I) 
+ (1 + (RE/4000 ) 2 )*RCMLS(2 , I) 
where RE = range to MLS elevation equipment 
TRI LATERAT I ON FILTER GAINS FLGTRI 
FLGTRI ( I ) * CONSTF( 10 + I) 

+ (1 + (R/500 ) 2 )*RCTRI ( 1 , I) 

+ (1 + R/1000)*RCTRI (2 , I) 
where R = total range from (0, 0, 0) 

ENROUTE FILTER GAINS FILGER 

FILGER(I) = CONSTF( 1 + I) + Z*RCER(1, I) 

+ (1 + V*Z*123.E-6)*RCER(2, I) 

where Z = estimated normal position 

(negative of altitude) 

V = estimated velocity magnitude 

FLADPT 
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TABLE 2 . 2-6a 

SPECIAL EVENT REPORT 


CMPUTR 


T-1724a 


REPORT 

DESCRIPTION - VARIABLES 

SUBROUTINE 

REPORT FROM 
SUBROUTINE 1 

CMPUTR 

TIME - Time of call to CMPUTR 

CMPUTR 


TABLE 2. 2 -6b 

SPECIAL EVENT REPORT - CNTRL1 


T-1724b 


' REPORT 

DESCRIPTION - VARIABLES 

SUBROUTINE 

REPORT FROM 

SUBROUTINE 

CNTRL1 

CNTRL1 FIRST CALL INITIALIZATION 

UVCCOM - Control command (true = estimate) 

YDTCOM - Guidance command push-down (control 
integrator state) 

CONTROL COMMAND IS - Control command to be 

sent to actuators 

CNTRL1 


TABLE 2.2-6C 

SPECIAL EVENT REPORT 


CNTRL2 


T-1724C 


REPORT 

DESCRIPTION - VARIABLES 

SUBROUTINE 

REPORT FROM 

SUBROUTINE 

CNTRL2 

CNTRL2 FIRST CALL 

UVCCOM - Control command (true = estimate) 
YDTCOM - Guidance command 

CONTROL COMMAND IS - Control command to be 

sent to actuators 

CNTRL2 


2-49 























THE ANALYTIC SCIENCES CORPORATION 


TABLE 2.2-7 

SPECIAL EVENT REPORT - GNADPT 


T-1725 


REPORT TITLE 

DESCRIPTION - VARIABLES 

SUBROUTINE 

REPORT FROM 
GNADPT 

Only issued if IGNPT = 1 (from NAMELIST/NINOU) 
Prints three gain scheduling variables, gains 

AIR RELATIVE VELOCITY MAGNITUDE- VT 

BODY Z-AXIS VELOCITY - W 

PSIDOT - <5 

STATE FEEDBACK MATRIX CX - 4x9 matrix 

LOW-PASS FILTER MATRIX CV - 4x4 matrix 

COMMAND INPUT MATRIX CY - 4x4 matrix 

COMMAND FEED- FORWARD MATRIX CYD - 4x4 matrix 

GNADPT 
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TABLE 2.2-8 

SPECIAL EVENT REPORT - FILINT 


T-1726 



REPORT TITLE 

DESCRIPTION - VARIABLES j 

SUBROUTINE 




REPORT FROM 
SUBROUTINE 

FILTER INITIALIZATION - Initial estimates 
set to true values (report only if KPRINT(2)>0) 

FILINT 


FILINT 

EARTH-RELATIVE POSITION - XIP/XIM 




EARTH-RELATIVE WIND - WIP/WIM 
BODY-AXES VELOCITY - VBP/VBM 
ACCELEROMETER BIAS - BBP/BBM 




BODY ANGULAR RATES - ARP /ARM 
EULER ANGLES - EAP/EAM 
AIR-RELATIVE VELOCITY - VAP/VAM 




EARTH-RELATIVE VELOCITY - VEP/VEM 
TRUE AIRSPEED - VTASP/VTASM 

SIDESLIP - BETAP/BETAM 


— 


MLS MEASUREMENTS - AMLSP/AMLSP 




IFLAGV - if 0 do not use air data 

if 1 use air data for all states 
if 2 use air data only to update wind 




IFLAGW - 0 propagate wind state 

1 do not propagate state 





K1 - if 1, transition between velocity 

> 50 kts and < 50 kts 




TIME1 - end of transition stage - < 50 kts to 
> 50 kts 


— 


TIME2 - end of transition stage - > 50 kts to 
< 50 kts 




VTAS1 - previous true airspeed estimate 
IFLGRM - MLS range flag 




HBE(3x3) - body to earth-axes transformation 
T 

HEB=(HBE) -earth-axes to body-axes transforma- 
tion 


— 


HME(3*3) - MLS-to-earth-axes transformation 
HEM(3*3) - (HME)” 1 earth-to-MLS transformation 




HBA(3x3) - body-to-air transformation 


— 


HAB(3x3) - (HBA) -1 air-to-body transfor- 
mation 




ALB( 3x3 ) - body angular rate to Euler angular 
rate transformation 


— 
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TABLE 2.2-9 

SPECIAL EVENT REPORT - FILTER 


T-1728 


REPORT TITLE 

DESCRIPTION - VARIABLES 

SUBROUTINE 

REPORT FROM SUB- 
ROUTINE FILTER 

If KPRINT(2 ) is set > 0 (by NAMELIST/NINOU / ) 
then full debug printout (Table 2.2-8) is 
produced during filter update and filter 
propagate. The following are printed in all 
cases : 

XEEST earth-relative position estimate 

(x, y, z) 

XVREST reduced state vector estimate 

(u, v, w, p, q, r, a, §, $) 

UVCEST estimated^ control 

(6 B , 6 c , 6 g , o R ) 

VEEST earth-relative velocity estimate 

AAA 

(X, y, z) 

VWEST earth-jelative wind estimate 

(wx , wy, wz) 

ZVCOBS observation vector from sensors 

ZVCOBS( 1-3 ) MLS or trilateration position 

ZVCOBS ( 4-6 ) trilateration velocity 

ZVCOBS (7-9 ) j accelerometer outputs ! 

ZVCOBS ( 10 ) barometric altimeter j 

ZVCOBS ( 11-13 ) Euler angle measurements 

ZVCOBS ( 14-16 ) body angular rates 

ZVCOBS (17) indicated airspeed 

ZVCOBS ( 18 ) sideslip 

ZVC0BS( 19-22 ) control position 

ZVCOBS (2 3-27 ) not used 

ZVCOBS(28-30 ) wind in body axes (true) 

ZVCOBS ( 31-42 ) true rigid-body state XVCTRU 
ZVC0BS( 31-46 ) true control command UVCTRU 

FILTER 

' 

1 

i 
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TABLE 2.2-10 

SPECIAL EVENT REPORT - GUIDE 


T-1729 


REPORT TITLE 

DESCRIPTION - VARIABLES 

SUBROUTINE 

REPORT FROM SUB- 
ROUTINE GUIDE 

GUIDANCE WAYPOINT NUMBER: 

INUM, the index into the prestored guidance 
trajectory such that current time, 

TIME > TIMG(INUM) 

TIMG - time at current guidance waypoint 

TIMG(INUM) 

TIME(NEXT)- time at next guidance waypoint 
TIMG (INUM + 1) 

DELTAT -TIMG(INUM + 1) - TIMG(INUM) : DT 

FAC: Fraction of current segment completed 

XIG - position (x, y, z) at waypoint INUM 

XI G( NEXT ) -position (x, y, z) at waypoint 
INUM + 1 

VIG - scheduled velocity (x, y, z) at 

waypoint INUM 

VIG (NEXT ) - scheduled velocity (x, y, z) at 
waypoint INUM + 1 

GUIDANCE GAINS - 4x4 guidance gain matrix at 
waypoint INUM 

EARTH-RELATIVE POSITION ERROR - DX difference 
between estimated position XEEST and scheduled 
guidance position XEG 

EARTH-RELATIVE VELOCITY ERROR - DV difference 
between estimated velocity VEEST and scheduled 
guidance velocity VEG 

EARTH-RELATIVE ACCELERATION ERROR - DVDOT 
difference between estimated velocity rate 1 

and guidance scheduled velocity rate J 

VELOCITY GUIDANCE COMMAND VECTOR - YDTCOM i 

(x c . y c , z c , j> c ) ; 

ATTITUDE GUIDANCE COMMAND VECTOR - YDTCOM i 

( 6 , <p , z . ip ) 1 

v c ' y c c ’ y c ' 

GUIDE 
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TABLE 2.2-11 

SPECIAL EVENT REPORT - COVPT 

T-1730 
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TABLE 2.2-11 (CONTINUED) 
SPECIAL EVENT REPORT - COVPT 


T-1731 


REPORT TITLE 

DESCRIPTION - VARIABLES 

SUBROUTINE 



CONTROL SDC 

X(I,I) = /P33(I,I) 


MATRIX 




Xfl.J) = P33 ^ 1 ' I f J 

COVPT 


y/P33( 1,1) P33(J,J) 

[ 


CONTROL ESTI- 

xti.j) = EMiLdl 

COVPT 

MATE CORRELA- 

yP33(I,I) P44 ( J , J ) 


TION MATRIX 



CONTROL CON- 

Xf T . J1 = PAW.J)', 

COVPT 

TROLLER CORRE- 

v/P33(I,I) P55( J , J) 


LATION MATRIX 



ESTIMATE SDC 

X( I , I ) = /P44(I,I) 


MATRIX 




Xfl.J) = F-A^L.J.L- I 4 j 

COVPT 


>44(1,1) P44( J , J) 


ESTIMATE CON- 

XfI r J) = 

COVPT 

TROLLER CORRE- 

>44(1.1) P45 ( J , J ) 


LATION MATRIX 



CONTROLLER SDC 

X(I,I) = >55(1.1) 

COVPT 

MATRIX 




X(I.J) = P55 ^' J) I 4 J 



>55(1,1) P55(J,J) 
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TABLE 2.2-12 
TIME HISTORY PR I NT /PLOT 


T-1732 


LA3EL 

DESCRIPTION 

H ' ■■ ■ ■■■ — — -■ ■ — 

SOURCE 

VARIABLE 

time 

Time of output sample 


TIME 

DEL-B 

Pitch control-differential collective 

6 b 

UVCTRU(l) 

DEL-C 

Vertical velocity control collective 

«c 

UVCTRU ( 2 ) 

XE 

Earth-relative axial position 

X 

UVCTRU ( 1 ) 

ZE 

Earth-relative normal position 

z 

XVCTRU ( 3 ) 

THETA 

Pitch Euler angle 

3 

XVCTRU(U) 

U 

Body x-axis velocity 

u 

XVCTRU ( 4 ) 

Q 

Body pitch angular rate 

q 

XVCTRU ( 8 ) 

w 

Body z-axis velocity 

w 

XVCTRU ( 6 ) 

V 

Body y-axis velocity 

V 

XVCTRU ( 5 ) 

R 

Body yaw angular rate 

r 

XVCTRU ( 9 ) 

P 

Body roll angular rate 

P 

XVCTRU ( 7 ) 

PHI 

Roll Euler angle 

$ 

XVCTRU (10) 

PSI 

Yaw Euler angle 

¥ 

XVCTRU (12) 

YE 

Earth-relative lateral position 

y 

XVCTRU ( 2 ) 

DEL-S 

Roll control-cyclic 

6 S 

UVCTRU ( 3 ) 

DEL-R 

Yaw control-differential cyclic 

6 r 

UVCTRU ( 4 ) 

XEDOT 

Earth-relative axial velocity 

X 

./ U \ 

YEDOT 

Earth-relative lateral velocity 

y 

M v 

ZEDOT 

Earth-relative normal velocity 

# 

Z 


DB-N 

Nominal pitch control-differential collective 

' 5 B 

UVNTRU ( 1 ) 

DC-N 

Nominal vertical velocity control-collective 

5 C 

UVNTRU ( 2 ) 

XE-N 

Nominal earth-relative axial position 

X 

XVNTRU ( 1 ) 

ZE-N ! 

Nominal earth-relative normal position 

Z 

XVNTRU ( 3 ) 

THETA-N 

Nominal pitch Euler angle 

e 

XVNTRU (11) 

U-N 

Nominal body x-axis velocity 

u 

XVNTRU (4) 

Q-N 

Nominal body pitch rate 

q 

XVNTRU ( 8 ) 

W-N 

Nominal body z-axis velocity 

w 

XVNTRU ( 6 ) 

i 

> 

Nominal body y-axis velocity 

V 

XVNTRU ( 5 ) 

R-N 

Nominal body yaw rate 

r 

XVNTRU ( 9 ) 
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TABLE 2.2-12 (CONTINUED) 
TIME HISTORY PRINT/PLOT 


T-X733 


LABEL 

DESCRIPTION 

SOURCE 

VARIABLE 


i 


P-N 

Nominal body roll rare 

P 

XVNTRU ( 7 ) 

PHI-N 

Nominal roll Euler angle 

0 

XVNTRU ( 10 ) 

PSI-N 

Nominal yaw Euler angle 

'V 

XVNTRU (12) 

YE-N 

Nominal earth-relative lateral position 

y 

XVNTRU ( 2 ) 

DS-N 

Nominal roll control-cyclic 

5 S 

UVNTRU ( 3 ) 

DR-N 

Nominal yaw control-differential cyclic 

Sr 

UVNTRU ( 4 ) 

DB-S 

Sensed control position 

°B 

ZVC0BS( 19 ) 

DC-S 

Sensed control position 

ic 

ZVCOBS ( 20 ) 

XE-EST 

Estimated earth-relative lateral position 

X 

XEEST(l) 

ZE-EST 

Estimated earth-relative normal position 

2 

XEEST ( 3 ) 

THET-EST 

Estimated pitch Euler angle 

9 

XVREST ( 8 ) 

U-EST 

Estimated body x-axis velocity 

u 

XVREST ( 1 ) 

Q-EST 

Estimated body pitch rate 

q 

XVREST ( 5 ) - 

W-EST 

Estimated body z-axis velocity 

w 

XVREST ( 3 ) 

V-EST 

Estimated body y-axis velocity 

V 

XVREST ( 2 ) 

R-EST 

Estimated body yaw rate 

r 

XVREST (6) 

P-EST 

Estimated body roll rate 

P 

XVREST (4) 

PHI-EST 

Estimated roll Euler angle 

A 

XVREST ( 7 ) 

PSI-EST 

Estimated yaw Euler angle 

i 

XVREST ( 9 ) 

YE-EST 

Estimated earth-relative lateral position 

y 

XEEST (2) 

DS-S 

| Sensed control position 

5 S 

ZVCOBS (21) 

DR-S 

Sensed control position 

5 R 

ZVCOBS (22) 

MLS-R 

"True" MLS range - uncorrupted by noise, - 

errors 

RAEXYZ 

MLS-AZ 

"True" MLS azimuth angle 


(XVCTRU( 1-3 ) ) 

MLS -EL 

"True" MLS elevation angle 



TRI-X 

"True" t rilaterat ion axial position 


XVCTRU(l) 

TRI-Y 

"True" t irlateration lateral position 


XV CTRL’ (2 ) 

TRI-Z 

"True" trilaterat ion normal position 


XVCTRU ( 3 ) 
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TABLE 2.2-12 (CONTINUED) 
TIME HISTORY PRINT/PLOT 


LABEL 

DESCRIPTION 

SOURCE 

VARIABLE 

VEXWND 

Actual earth-relative axial wind velocity 

AVINDB(l) 

VEYV.'ND 

Actual earth-relative lateral wind velocity 

Hg| WINDB( 2 ) 
'WINDB( 3 ) 

VEZWND 

Actual earth-relative normal wind velocity 

BIAS-U 

Accelerometer bias - body x-axis b u 

SVCTRU ( 5 ) 

BIAS-V 

Accelerometer bias - body y-axis b v 

SVCTRU ( 6 ) 

BIAS-W 

Accelerometer bias - body z-axis b w 

SVCTRU ( 7 ) 

DBCOM 

Pitch control command from controller 

UVCCOM(l) 

DCCOM 

Vertical velocity command from controller 

UVCCOM( 2 ) 

DSCOM 

Roll control from controller 

UVCC0M(3) 

DRCOM 

Yaw control command from controller 

UVCCOM(4) 

YD1 

Controller integrator state (guidance command 1) 

YDTCOM ( 1 ) 

YD2 

Controller integrator state (guidance command 2) 

YDTC0M(2) 

YD 3 

Controller integrator state (guidance command 3) 

YDTCOM ( 3 ) 

YD4 

Controller integrator state (guidance command 4) 

YDT COM ( 4 ) 

DBACT 

Actuator pitch control dynamic state o BA 

XVCTRU (21) 

DCACT 

Actuator vertical velocity control dynamic 5 a , 

state u 

! XVCTRU (22) 

DSACT 

Actuator roll control dynamic state o SA 

XVCTRU (23) 

DRACT | 

Actuator yaw control dynamic state 

XVCTRU (24) 

DBDOT 

Differential collective rate pitch control rateoA ! 
rate ' 

XVCTRU (17) 

DCDOT 

Collective rate - vertical velocity control oa 

rate 

XVCTRU (IS) 

DSDOT 

Cyclic rate - roll control rate $2 

XVCTRU (19) 

DRDOT 

Differential cyclic rate - yaw control rate 

XVCTRU (20) 
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TABLE 2.2-12 (CONTINUED) 
TIME HISTORY PRINT/PLOT 


T-173o 

LABEL 

DESCRIPTION 


SOURCE 

VARIABLE 

ACC-XB 

Accelerometer output -body y-axis 

AU 

ZVC0BS( 7 ) 

ACC-YB 

Accelerometer output -body y-axis 

AV 

ZVC0BS(8 ) 

ACC- ZB 

Accelerometer output-body z-axis 

AW 

ZVC0BS(9 ) 

ALT 

Barometric altimeter output 

Z B 

ZVCOBS(IO) 

ROLL-S 

Sensed roll Euler angle 

* 

ZVCOBS(ll) 

PITCH-S 

Sensed pitch Euler angle 

e 

ZVC0BS( 12 ) 

YAW-S 

Sensed yaw Euler angle 

iy 

ZVC0BS( 13 ) 

P-S 

Sensed body roll rate 

P 

ZVC0BS( 14 ) 

Q-S 

Sensed body pitch rate 

■ q 

ZVC0BS(15 ) 

R-S 

Sensed body yaw rate 

r 

ZVC0BSC16) 

IAS 

Indicated airspeed 

IAS 

ZVC0BS( 17 ) 

BETA 

Sensed sideslip 

3 

ZVC0BS( 18 ) 

XEDT-EST 

Estimated earth-relative axial velocity 

X 

VEEST(l) 

YEDT-EST 

Estimated earth-relative lateral velocity 

y 

VEEST( 2 ) 

ZEDT-EST 

Estimated earth-relative normal velocity 

z 

VEEST( 3 ) 

XEDOT-G 

Scheduled guidance axial velocity 

*d 

VEG(l) 

VEDOT-G 

Scheduled guidance lateral velocity 

^G 

VEG(2 ) 

2EDOT-G 

Scheduled guidance normal velocity 

*0 

VEG( 3 ) 
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TABLE 2.2-12 (CONTINUED) 
TIME HISTORY PRINT /PLOT 


T-1736 


LABEL 

DESCRIPTION 


SOURCE 

VARIABLE 

XE-G 

Scheduled guidance axial position 

X G 

XEG(l) 

YE-G 

Scheduled guidance lateral position 

y G 

XEG(2 ) 

ZE-G 

Scheduled guidance normal position 

Z G 

XEG( 3 ) 

K(Q-DR) 

Control feedback gain - pitch rate to yaw control 

CX(4, 5) 

K(R-DS) 

Control feedback gain - yaw rate to roll control 

CX( 3 , 6) 

S VX-DB 

Control input gain - axial velocity to pitch 
. control 


CY(1, 1) 

D VY-DS 

Control feedforward gain - lateral velocity to 
roll control 


CYD ( 3 , 2) 

RNG-EST 

Estimated range to MLS equipment 

R 

XP( 27 ) 

AZ-EST 

Estimated azimuth to MLS azimuth tower 

a 

XP ( 2 8 ) 

EL-EST 

Estimated elevation angle to MLS elevation 
tower 

A 

e 

XP( 29 ) 

VTAS-EST 

Estimated air speed 

V a 

XP( 25 ) 

BETA-EST 

Estimated sideslip 

s 

XP ( 2 6 ) 

WEX-EST 

Estimated axial wind-earth-relative 


VWEST(l) 

WEY-EST 

Estimated lateral wind-earth-relative 

?y 

w z 

VWEST( 2 ) 

WEZ-EST ; 

Estimated normal wind-earth-relative 

VV/EST ( 3 ) 

B-U-EST ' 

Estimated accelerometer bias - body x-axis 

a" 

b U 

XP ( 7 ) 

B-V-EST 

Estimated accelerometer bias - body y-axis 

b v 

XP ( 8 ) 

B-W-EST 

VTAS 

Estimated accelerometer bias - body z-axis 
True airspeed 

A ¥ | 

b w 

XP ( 9 ) 

BETA-T 

True sideslip 


i 
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TABLE 2.2-12 (CONTINUED) 
TIME HISTORY PRINT /PLOT 



LABEL 


DESCRIPTION 

RMS : 

XE 

Standard deviation 

of helicopter axial position 

RMS : 

YE 

Standard deviation 

of helicopter lateral position 

RMS: 

2E 

Standard deviation 

of helicopter normal position 

RMS : 

U 

Standard deviation 

- body x-axis velocity 

RMS: 

V 

Standard deviation 

- body y-axis velocity 

RMS: 

w 

Standard deviation 

- body z-axis velocity 

RMS: 

p 

Standard deviation 

- body roll rate 

RMS: 

Q 

Standard deviation 

- body pitch rate 

RMS : 

R 

Standard deviation 

- body yaw rate 

RMS: 

PHI 

Standard deviation 

- roll Euler angle 

RMS : 

THET 

Standard deviation 

- pitch Euler angle 

RMS : 

PSI 

Standard deviation 

- yaw Euler angle 

RMS : 

uw 

Standard deviation 

- turbulence on body x-axis 

RMS: 

vw 

Standard deviation 

- turbulence on body y-axis 

RMS : 

ww 

Standard deviation 

- turbulence on body z-axis 

COR: 

u-uw 

Correlation - body 
turbulence 

x-axis velocity 'to body x-axis 

COR: 

v-vw 

Correlation - body 
turbulence 

y-axis velocity to body y-axis 

COR: 

w-ww 

Correlation - body 
turbulence 

z-axis velocity to body z-axis 

COR: 

P-VW 

Correlation - body 
turbulence 

roll rate with body y-axis 

COR: 

Q-WV 

Correlation - body 
turbulence 

pitch rate with body z-axis 

COR: 

R-VW 

Correlation - body 
turbulence 

yaw rate with body y-axis 

CCR : 

TH-U 

Correlation - body 
x-axis velocity 

pitch Euler angle with body 

COR: 

W-Q 

Correlation - body 
pitch rate 

z-axis velocity with body 

COR: 

u-w 

Correlation - body 
z-axis velocity 

x-axis velocity with body 

COR: 

V-P 

Correlation - body 

y-axis velocity with roll rate 

COR: 

V-R 

Correlation - body 

y-axis velocity with yaw rate 

COR: 

v-w 

Correlation - body 
z-axis velocity 

y-axis velocity with body 


SOURCE 

VARIABLE 


v'PIKI.l) 


%/PH(2,2) 


y/Pll(3,3) 


v/Pll(4,4) 


v'Pll(5,5) 


Pll( 6,6) 


v/Pll (7 , 7 ) 


Pll(8,8) 



v/Pll(10,10) 


P 11(11, 11) 



v/Pll(6,6) P22 (3,3) 
P12(7 . 2 ) 


,/Pll(7,7) P22(2,2) 
P12C8.3) 


v/Pll (8,8) P22 (3,3) 
P12 (9,21 


\/Pll( 9 , 9 ) P22(2 ,2) 
P1K4.11) 



v/Pll(4,4) Pll(6,6) 
Pll (5,7 
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TABLE 2.2-12 (CONTINUED) 
TIME HISTORY PRINT/PLOT 


T-1738 


1 

LABEL | 

DESCRIPTION 

SOURCE 

VARIABLE 

. ; - - . . 1 


COR: 

Q-R 

Correlation - pitch rate with yaw rate 

Pll( 8 , 9 ) 

v/Pll( 8 , 8 ) Pll(9,9) 

COR: 

Q-DB 

Correlation - pitch rate with pitch control 

P13( 8 , 1 ) 

v/?ll(8,8) P33( 1 , 1 ) 

COR: 

W-DC 

Correlation - body z-axis velocity with vertical 
control 

P13 ( 6,2) 

>11(6,6) P33(2 , 2 ) 

COR: 

V-DS 

Correlation - body y-axis velocity with roll 
control 

P13(5 . 3) 

>11(5,5) P33( 3 , 3 ) 

COR: 

P-DS 

Correlation - roll rate with roll control 

P13(7.3) 

>11(7 ,7) P33( 3 , 3) 

COR: 

R-DR 

Correlation - yaw rate with yaw control 

P13( 9 , 4 ) 

n/P11(9,9) P33(4 , 4 ) 

RMS 

DB 

Standard deviation - pitch control 

>33(1.1) 

RMS 

DC 

Standard deviation - vertical velocity control 

/P33(2 ,2 ) 

RMS 

DS 

Standard deviation - roll control 

>33(3,3) 

RMS 

DR 

Standard deviation - yaw control 

>33(4,4) 

RMS 

XEST 

Standard deviation - axial position estimate 

>44 (1,1) 

RMS 

YEST 

Standard deviation - lateral position estimate 

>44 (2,2) 

RMS 

ZEST 

Standard deviation - normal position estimate 

>44 (3,3) 

RMS 

UEST 

Standard deviation - body x-axis velocity estimate 

>44 (4,4) 

RMS 

VEST 

Standard deviation - body y-axis velocity estimate 

>44(5,5) 

RMS 

TOST 

Standard deviation - body z-axis velocity estimate 

>44(6,6) 

RMS 

PEST 

Standard deviation - estimated roll rate 

>44 (7,7) 

RMS 

QEST 

Standard deviation - estimated pitch rate 

/P44 (8,8) 

RMS 

REST 

Standard deviation - estimated yaw rate 

> 44 ( 9 , 9 ) 

RMS 

PHIE 

Standard deviation - estimated roll Euler angle 

>44(10,10) 

RMS 

THEE 

i 

Standard deviation - estimated pitch Euler angle 

>44(11,11) 

RMS 

PSIE 

Standard deviation - estimated yaw Euler angle 

>44(12,12) 

RMS 

INTI 

1 

Standard deviation - guidance commanded axial 
velocity 

>55(1,1) 

RMS: 

: INT2 | 

Standard deviation - guidance commanded lateral 
velocity 

>55(2,2) 

RMS: 

: INT3 

i 

i 

i 

Standard deviation - guidance commanded normal 
velocity 

>55(3,3) 

RMS' 

: INT4 

Standard deviation - guidance commanded yaw Euler 
angle 

>55(4,4) 

P35Q.1) 

>33(1,1) P55(l ,1) 

CO: 

DB- 1 1 

Correlation - pitch control with axial velocity 
command 

CO: 

DC- 11 

Correlation - vertical velocity control with 
lateral velocity command 

P35 ( 2.2) 

>33(2,2) P55 (2,2) 

CO: 

IDS- I 3 

Correlation - roll control with normal velocity 
command 

P35( 3.3) 

>35(3,3) P55( 3,3) 

CO. 

DIR-I4 : 

Correlation - yaw control with yaw angle command 

P35 ( 4.4) 


i — : 

i/P33(4,4) P55 ( 4 . 4 ) 
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APPENDIX A 
UTILITY ROUTINES 

The following utility routines are used by TVHIS 
to accomplish the described functions. 

A . 1 MATRIX UTILITY ROUTINES 

A. 1.1 ADD Computes the sum of 2 matrices. 

A. 1.2 MNUS Sets a matrix to the negative of a given 

matrix. 

A. 1.3 MULT Computes the product of 2 matrices. 

A. 1.4 INV Inverts a general matrix by Gauss- Jordan 

elimination. 

A. 1.5 SUBT Computes the difference of 2 matrices. 

A. 1.6 TRNP Computes the transpose of a matrix. 

A. 1.7 EQAT Sets a matrix equal to another matrix. 

A. 1.8 SMLT Multiplies a matrix by a scalar. 

A. 1.9 I DEN Sets a square matrix to the identity matrix. 


A-l 



THE ANALYTIC SCIENCES CORPORATION 


A. 1.10 ZERO Sets all elements of a matrix to zero. 

A. 1.11 CPEQ Computes the cross-product-equivalent 

matrix for a vector. 

A. 1.12 NORMAL Normalizes a given covariance matrix. 

A. 1.13 RESTOR Recovers a symmetric matrix from its 

lower-triangular representation in a 
vector. 

A. 1.14 STORE Stores the lower-triangular representation 

of a symmetric matrix in a vector. 

A. 1.15 TAPP Computes the sum of a matrix with the 

transpose of a second matrix. 

A. 1.16 TMLT Computes the product of a matrix with 

the transpose of a second matrix. 

A. 2 COORPINATE TRANSFORMATION ROUTINES 

A. 2.1 HBAS Calculates body axes to air-relative axes 

Transformation matrix. 

A. 2. 2 HBES Calculates body axes to earth axes 

transformation matrix. 

A. 2. 3 HBIR Calculates general body axes to inertial axes 

Transformation matrix. 

A. 2. 4 HEMS Calculates earth axes to MLS axes trans- 
formation matrix. 
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A. 2. 5 LB I Calculates body angular rates to Euler 

angular rates transformation matrix. 

A. 2. 6 LBPR Calculates partial derivative of body angular 

rates with respect to Euler angle. 

A. 2. 7 LBR Calculates Euler angular rates to body 

angular rates transformation matrix. 

A. 2. 8 RAEXYZ Transforms earth-relative position to 

MLS range, azimuth and elevation. 

A. 3 UTILITY PRINT ROUTINES 


A. 3.1 

PRMR 

Prints a 

real-valued matrix in 

tabular form. 

A. 3.2 

WFHPLT 

Produces 

printer-plots . 


A. 3. 3 

COVPT 

Prints a 

normalized covariance 

matrix . 


A. 4 RANDOM-NUMBER GENERATION AND TRIGONOMETRIC FUNCTIONS 

A. 4.1 GAUSS Generates a normally distributed pseudo- 
random number. 

A. 4. 2 RANDOM Generates a uniformly distributed pseudo- 
random number. 

A. 4. 3 TRIGS Calculates second-order Taylor series 

approximations for sine, cosine, tangent, 
and secant functions. 
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APPENDIX B 


IMPLEMENTATION OF PERMANENT FILE FOR 
AERODYNAMIC COEFFICIENTS 


The program (AERO) must be run before TVHIS to store 
aerodynamic coefficients on a permanent disk file. The file 
consists of 297 different cases of aerodynamic coefficients 
whose parameter variations are explained below. 

The aerodynamic coefficients, corresponding to specific 
values of ip, x and z are selected by computing a case number 
(NCASE) defined as: 

NCASE = ( J-l)*99 + (K-l )*11 + L 

J is an index of 3 values for ip, K is an index of 9 values for 
z, and L in an index of 11 values for x. The ij) (radians/sec) 
values of 0.0, +0.05, -0.05 correspond to J = 1 to 3. The z 
(feet/min) values of -2000, -1500, -1000, -500, 0, 500, 1000, 
1500, 2000 correspond to K = 1 to 9. The x (knots) values 
of -40, -20, 0, 20, 40, 60, 80, 100, 120, 140, 160 correspond 
to L = 1 to 11. For each case, there are 71 aerodynamic co- 
efficients used in the formulation of the system matrix F. 
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APPENDIX C 

LIST OF DELIVERABLES 


Delivered with this document are nine distinct, 
computer-related items: one magnetic tape, four executable 

card decks, and four computer-printed listings. These nine 
items are detailed in the following list. 


1. Magnetic tape, 7-track, 556 bpi (HI) con- 
taining three files: 

File 1 (binary) TVHIS program library 
File 2 (binary) AERO program library 
File 3 (coded) aerodynamic coefficient data 

2. Executable deck which reads Item 1 (above) 
and creates two permanent files. 

TVHOBJ - TVHIS object code 
AEROB - data file of aerodynamic coeffi- 
cients, read by TVHIS 

3. Executable deck - TVHIS sample execution # 1 

4. Executable deck - TVHIS sample execution #2 

5. Executable deck - TVHIS sample execution #3 

6. Computer printout generated by execution 

of Item 2 — includes compiler-generated 
source listings of TVHIS and AERO, and the 
execution time output of program AERO 


7. 


Computer printout generated by execution 
of Item 3 - TVHIS sample output #1 


8 . 


Computer printout generated by execution 
of Item 4 - TVHIS sample output #2 


9. 


Computer printout generated by execution 
of Item 5 - TVHIS sample output #3 
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